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Abstract 

This  investigation  of  rotor-stator  interaction  noise  originated  from  an  interest  in 
reducing  a  ventilation-noise  problem.  This  type  of  noise  is  generated  aerodynamically  by 
the  flow  interaction  between  the  rotor  and  stator.  In  order  to  reduce  this  interaction  noise 
thoroughly  along  the  pipe,  an  efficient  method  is  to  remove  noise  sources 
aerodynamically  in  the  interaction  flow  region  by  means  of  active  flow  control.  This 
method  suggests  that  we  need  to  identify  the  active  noise  source  region  before  attempting 
to  implement  the  noise  reduction.  On  this  account,  we  shall  employ  measurement 
techniques  for  this  identification  process  by  correlating  both  the  interaction-flow  field 
(cause)  and  the  corresponding  acoustic  field  (effect).  Therefore,  the  motive  of  the  present 
study  is  to  provide  a  general  basis  for  locating  the  noise  sources  experimentally  and  to 
pave  the  road  for  possible  application  of  the  active  rotor-stator  interaction-noise  control 
in  the  future. 

1.1.  Theories  regarding  the  Present  Study 

The  theory  of  'aeroacoustics’  or  'sound  generated  aerodynamically’  is  built  upon 
the  governing  equations  of  a  fluid  flow,  i.e.,  the  continuity  and  momentum  equations. 
There  are  two  approaches  to  the  study  of  aeroacoustics  (Goldstein  1976):  one  is  to 
directly  solve  the  linearized  governing  equations  specifically  for  sound  generation  due  to 
pressure  fluctuations  on  a  solid  boundary  in  a  moving  medium  (Kaji  and  Okazaki  1970, 
Myers  and  Kershen  1995),  the  other  is  to  follow  Lighthill’s  acoustic  analogy  (Lighthill 
1952,  1954,  1962)  deduced  from  the  governing  equations  (Powell  1960,  Meecham  1965, 
Ffowcs  Williams  and  Hawkings  1969a,  1970,  Fujita  and  Kovasznay  1974).  We  will 
follow  the  second  approach  to  investigate  the  noise  generated  by  rotor-stator  interaction. 

Based  on  Lighthill’s  theory,  Ffowcs  Williams  and  Hawkings  (1969b)  have 
established  a  more  general  aeroacoustic  theory  for  a  moving  solid  boundary  in  a 
nonstationary  medium.  Their  theory  provides  a  basis  for  noise  generation  due  to  rotor- 
stator  interaction.  However,  in  the  reported  experimental  papers,  a  transient  stator 
response  is  usually  responsible  for  the  majority  of  the  interaction-noise  generation  (Fujita 
and  Kovasznay  1974,  Ho  and  Kovasznay  1976a,  Zandbergen  1988,  Simonich  et  al. 
1993).  Since  the  stator  is  stationary  in  the  flow,  we  may  simplify  Ffowcs  Williams- 
Hawkings  theory  in  order  to  obtain  a  relationship  between  the  stator  response  and  the 
noise  generation.  The  sound  field  radiated  by  a  stationary  solid  body  in  a  moving 
medium  (Fig.  1-1)  is 
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where  i?  =  |x-y|.  Eq.(l.l)  was  first  derived  by  Curie  (1955)  based  on  Lighthill’s 
theory. 

The  volume  integral  in  Eq.(l.l)  is  exactly  the  solution  of  Lighthill’s  equation 
(Lighthill  1952)  and  represents  a  sound  field  produced  by  a  volume  distribution  of 
quadrupoles.  Similarly,  the  surface  integral  depicts  an  equivalent  system  to  the  sound 
generated  in  a  stationary  medium  by  a  distribution  of  dipoles  of  strength  f.  per  unit  area, 
where/;  is  the  fluctuating  force  per  unit  area  exerted  by  the  solid  boundaries  on  the  fluid 
in  the  x.  direction.  Therefore,  the  surface  integral  term  is  known  as  Curie’s 
supplementary  formula  to  the  sound  field  with  the  existence  of  a  stationary  solid 
boundary  in  a  moving  medium  (Clark  and  Ribner  1969,  Siddon  1973,  Fujita  and 
Kovasznay  1974). 

It  should  be  noticed  that  t  -  R/a^,  the  so-called  'retarded  time  ’,  represents  the  time 
at  which  the  sound  received  at  time  r  by  a  sensor  at  the  point  x  was  radiated  from  the 
source  at  the  point  y,  where  /?  =  lx  -  yl,  as  illustrated  in  Fig.  1-1.  In  other  words,  an  event 
(sound)  happening  at  time  t  -  shall  be  observed  at  time  t  by  an  observer  (sensor) 
located  at  a  distance  R. 

By  neglecting  the  quadrupole  term  based  on  an  order-of-magnitude  comparison. 
Curie  (1955)  imposed  a  far-field  condition  on  Fq.(l.l)  and  came  to  an  estimate  of  the 
far-field  noise  or  radiation  field  for  low  Mach-number  flow  with  the  presence  of  a  solid 
body  in  the  form 

(1.2) 
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unit  area  exerted  by  the  solid  surface  on  the  fluid. 

Therefore,  Fq.(1.2)  is  a  simplified  Curie’s  supplementary  formula  for  the  far-field  noise. 

Clark  and  Ribner  (1969)  have  proven  the  validity  of  Fq.(1.2)  experimentally  by 
finding  the  direct  correlation  between  far-field  sound  and  fluctuating  lift  on  an  airfoil  in 
turbulent  flow.  In  parallel  with  the  work  of  Clark  and  Ribner,  Siddon  (1973)  established 
a  formal  cross-correlation  method,  from  Fq.(1.2),  to  calculate  surface  dipole  strength 
from  both  the  unsteady  pressure  on  an  airfoil  in  a  jet  flow  and  the  corresponding  far-field 
noise.  Siddon ’s  result  indicates  that  the  distribution  of  surface  dipole  strength  is 
identified  locally  on  the  solid  boundary  and  does  not  necessarily  correspond  to  the 
pressure  distribution  on  the  same  boundary.  This  is  because,  from  Fq.(1.2),  only  the  dFf 
dt  part,  i.e.,  the  time  derivative  of  pressure  fluctuations,  not  F.  itself,  contributes  to  the 
far-field  noise.  However,  for  those  measuring  positions  in  the  near  field,  Fq.(1.2)  is  no 
longer  satisfied.  The  near-field  noise  should  be  obtained  by  reconsidering  Fq.(l.l). 
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For  the  near  field,  Fujita  and  Kovasznay  (1974)  defined  a  less  restrictive 
condition,  that  is,  only  the  square  of  the  sound-propagating  distance  is  much  larger  than 
that  of  the  dimensions  of  the  source  region.  Under  this  assumption,  they  were  able  to  re¬ 
derive,  from  Curie’s  supplementary  formula,  a  simplified  equation  for  noise  prediction  in 
the  near  field.  They  found  that  the  results  obtained  from  noise  prediction  were  in  good 
agreement  with  those  measured  from  experiments.  This  simplified  equation  for  free- 
space,  near-field  dipole  sound  is 
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where 


is  the  freestream  dynamic  pressure  with  freestream  speed  , 


0050  =  Xj  /  X  the  directivity  of  the  dipole  sound,  R  /  a^  the  time  delay  due  to  the 
propagating  distance  of  sound  wave  between  source  and  detector,  A  the  area  of  stator, 
and  Q  the  lift  coefficient  of  the  stationary  body.  In  fact,  this  simplification  for  the  near- 
field  noise  has  been  already  examined  by  Lighthill  (1962)  in  an  explanation  of  sound 
generation  by  a  simple  dipole  in  a  non-moving  medium.  Similar  to  Lighthill’ s 
description,  the  first  term  on  the  right-hand  side  of  Eq.(1.3),  with  Q  which  represents  the 
total  lift  fluctuation,  is  identified  as  the  near-field  noise,  while  the  second  term,  the 
radiation  field  or  far-field  noise,  of  the  same  as  that  shown  in  Eq.(1.2). 

In  light  of  the  above  theoretical  results,  we  will  be  able  to  investigate  the 
interaction  noise  from  a  rotor-stator  stage.  Before  delving  deeper  into  the  current 
research  interest,  a  brief  discussion  will  be  made  about  the  features  of  the  rotor-stator 
interaction  noise. 


1.2.  Characteristics  of  Rotor-Stator  Interaction  Noise 


A  typical  sound-pressure-level  spectrum  of  the  microphone  signal  of  rotor-stator 
interaction  noise  is  shown  in  Fig.  1-2,  for  an  arrangement  with  a  two-blade  rotor  and  a 
two-blade  stator.  The  blade  passing  frequency  (BPF)  is  50Hz  with  a  rotor  rotation  of 
25Hz,  and  is  identified  by  a  clear  spike  at  50Hz  in  Fig.  1-2.  We  also  observe,  from  Fig.l- 
2,  that  there  exist  discrete  spikes  at  100,  150,  200  and  250  Hz,  respectively,  which  shows 
that  the  sound  radiation  due  to  the  rotor-stator  interaction  is  at  discrete  frequencies  that 
are  multiples  of  the  blade  passing  frequency  (Kaji  and  Okazaki  1970,  Hanson  1973, 
Goldstein  1976,  Blake  1986). 

Since  the  wave  energy,  in  Fig.  1-2,  is  mainly  distributed  at  the  multiple  BPF’s, 
the  signal  in  the  time  domain  should  be,  by  definition,  complex  periodic  data  whose 
waveform  repeats  itself  at  regular  intervals  (Bendat  and  Piersol  1986).  In  addition,  the 
signals  measured  from  the  flow  properties  (for  example,  surface  pressure,  velocity  field, 
etc.)  around  the  rotor-stator  stage  are  also  complex  periodic  due  to  the  rotor’s  rotation.  A 
phase-averaging  technique,  hence,  is  applicable  to  the  current  study,  and  will  be 
introduced  in  Chapter  2. 
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After  the  sound  radiates  from  the  rotor-stator  interaction  region,  it  propagates 
through  the  pipe  in  both  directions,  upstream  and  downstream.  For  the  sound  propagating 
downstream,  a  vortical  shear  flow  induced  by  the  rotor’s  rotation  will  refract  the  sound 
from  its  propagating  path  and  may  also  attenuate  it  a  little  by  heat  or  viscous  dissipation 
(Ho  and  Kovasznay  1976b,  Goldstein  1976).  The  theory  of  ray  tracing  for  geometric 
acoustics  is  devoted  to  such  a  wave  propagation  through  a  shear  flow  (Pridmore-Brown 
1962,  Thompson  1971,  Lighthill  1978,  Pierce  1989).  When  flow  velocities  are  much 
smaller  than  the  sound  speed,  the  interaction  between  the  sound  wave  and  the  vortical 
shear  flow  is  involved  with  very  gradual  changes  (Lighthill  1978).  Since  in  this  study  the 
flow  is  incompressible,  this  interaction  hence  will  be  ignored. 

Another  point  of  special  interest  in  the  rotor-stator  interaction  noise  is  regarding 
the  noise  resource  in  the  interaction-flow  field.  The  transient  stator  response  due  to  the 
passing  wake  of  the  rotor  used  to  be  deemed  as  the  noise  resource  (Fujita  and  Kovasznay 
1974,  Ho  and  Kovasznay  1976a,  Zandbergen  1988,  Simonich  et  al.  1993).  Nevertheless, 
still  very  little  is  known  about  the  details  of  how  this  noise  source  is  influenced  by  the 
interaction-flow  field.  Therefore,  the  current  investigation  will  merely  employ  the 
experimental  methods  to  identify  the  noise-source  region  from  the  associated  interaction- 
flow  field. 

The  gap  between  the  rotor  and  stator  also  plays  a  very  important  role  on  the 
interaction-noise  generation  (Sugeng  and  Fiedler  1986).  We  may  expect  that  the  shorter 
the  gap,  the  stronger  the  flow  interaction,  and  hence  the  louder  the  noise.  For  simplicity, 
we  shall  confine  a  fixed  gap  for  investigating  the  noise-generation  mechanism  in  this 
study. 

Inasmuch  as  this  research  interest  stems  from  analyzing  the  noise  generation  in  a 
ventilation  system,  the  flow  concerned  in  the  current  study  will  be  incompressible. 
Nevertheless,  the  results  should  be  helpful  as  a  basis  for  elucidating  the  mechanism  of  the 
interaction-noise  generation  in  a  more  complicated  turbomachinery  system. 
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1.3.  Goal  of  the  Present  Study 


The  goal  of  this  research  work  is  to  understand,  by  experiments,  the  physical 
mechanism  of  noise  generation  due  to  rotor-stator  interaction  in  an  incompressible  pipe 
flow.  The  approach  of  this  experimental  study  is  attempting  to  correlate  the  measured 
acoustic  field  at  a  certain  cross-section  of  the  pipe  to  the  aerodynamic  field  observed  in 
the  flow  region  of  the  rotor-stator  interaction.  There  are  three  different  measurements 
involved  with  this  investigation;  they  are  downstream  noise  measurements,  surface 
pressure  measurements  on  the  stator,  and  velocity  measurements  around  the  stator.  Since 
the  raw  data  from  all  measurements  are  periodic  based  on  BPF,  the  results  from  these 
data  will  be  finalized  by  a  phase-averaging  technique. 

After  this  experimental  investigation,  we  should  be  able  to  determine  the  active 
interaction  phase  interval,  to  identify  the  active  noise  source  region  in  the  interaction- 
flow  field,  and  to  understand  the  physical  mechanism  of  noise  generation  from  the  rotor- 
stator  interaction. 
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CHAPTER  2 

EXPERIMENTAL  CONFIGURATION 


2.1.  Wind  Tunnel  and  Test  Section 

A  blow-down  open-loop  wind  tunnel  with  an  anechoic  settling  chamber  was 
employed  for  the  experimental  setup  (Fig.  2-1).  The  maximum  freestream  speed  is  about 
15  m/s.  The  test  section  is  a  circular  plexiglass  pipe  with  dimensions  of  0.7  m  in  length 
and  a  0.178  ra  inner  diameter  (7"  I.D.).  A  detachable  plexiglass  section  is  connected  to 
the  rear  end  of  the  test  section  for  microphone  measurements,  as  illustrated  in  Fig.  2-1. 

Rotor  and  Stator 

For  simplicity,  a  test  stage  consisting  of  a  two-blade-rotor  and  a  two-blade-stator 
combination  in  a  circular  pipe  has  been  set  up  for  this  study  (Fig. 2-1  &  2-2).  The  stator 
was  located  at  a  distance  0.71  m  upstream  as  measured  from  the  end  of  the  rear  pipe 
section  to  the  leading  edge  of  the  stator  blade  (Fig.2-1). 

Instead  of  a  circular  cylinder  for  the  rotor  which  was  commonly  used  in  the 
reported  experimental  studies  (Fujita  and  Kovasznay  1974,  Ho  and  Kovasznay  1976a, 
Sugeng  and  Fiedler  1986,  Zandbergen  1988,  Simonich  et  al  1993),  a  fan-type-profile 
propeller,  with  variation  in  both  chord-length  and  pitch  along  the  span,  has  been  adopted 
as  the  rotor,  which  facilitates  a  variety  of  operating  conditions  for  investigation. 

Accordingly,  a  6-inch-pitch,  7-inch-diameter,  2-blade  model-plane  propeller 
(from  the  Master  Airscrew  Series  made  by  Windsor  Propeller  Co.)  is  chosen  for  the 
rotor.  The  rotor  was  driven  by  a  DC  motor  and  operated  below  4000  RPM.  In  contrast  to 
the  selection  of  the  rotor,  the  stator  consists  of  two  constant-NACA0016-profile  blades 
with  a  small,  bronze  hub  (with  0.026m  O.D.  and  0.035m  in  length)  between  them.  Each 
stator  blade  has  a  chord  length  of  0.05  m  and  a  span  of  0.076  m.  Comparing  the 
dimensions  of  the  stator  blades  and  the  hub,  it  can  be  asserted  that  the  noise  generated  by 
the  hub  should  be  much  smaller  than  the  noise  from  the  rotor-stator  interaction  in  this 
experimental  setup.  Hence,  the  noise  generated  by  the  rotor-stator  interaction  has  been 
mainly  investigated  in  this  study.  The  pressure  side  of  the  stator  is  defined  as  the  side 
facing  the  approaching  rotor,  and,  in  turn,  the  suction  side  of  the  stator  is  on  the  leeward 
side  (Fig.2-3). 

Definitions  of  phase  angle,  at  0°and  ±180° 

The  output  of  a  micro-switch  sensor  attached  on  the  rotor's  driving  motor  is  used 
as  a  phase  reference  signal  based  on  the  rotor's  rotation.  With  this  reference,  a  phase¬ 
averaging  technique  can  be  applied  to  process  raw  data  from  three  different 
measurements,  and  will  be  described  in  Section  2.3.1. 

The  reference  signal  is  a  pulse  related  to  the  instant  at  which  the  rotor's  trailing 
edge  passes  in  front  of  the  stator's  leading  edge.  Since  the  stator's  angle  of  attack  (AoA) 


6 


is  changed  from  case  to  case  in  this  study,  this  signal,  therefore,  is  triggered  relative  to 
the  stator's  position  at  zero-degree  AoA,  as  shown  in  Fig.2-2(b).  With  this  arrangement, 
we  will  be  able  to  define  a  time  reference  named  zero  degree  phase  angle,  i.e.,  </>  =  0°,  in 
order  to  assist  data  analysis.  The  phase  angle  (j>  at  4-/- 180°  is  defined  as  the  rotor  rotating 
to  a  vertical  position  relative  to  the  stator,  as  illustrated  in  Fig.  2-2(b). 

Stator’s  angle  of  attack,  a, ,  &  rotor’s  effective  angle  of  attack  ,  fie 

With  respect  to  the  freestream  direction,  the  stator's  angle  of  attack  (AoA ),  a,,  is 
in  a  positive  sense  when  the  stator  is  turned  to  its  pressure  side,  in  Fig.  2-3. 

An  effective  angle  of  attack  ,  fie,  of  a  rotor-blade  section  at  a  certain  radius  of  the 
rotor  is  determined  by  the  freestream  velocity  ,  the  rotor's  velocity  and  the  zero- 
lift-line  angle  Oof  the  same  rotor  section,  as  illustrated  in  Fig.2-3.  This  relationship  is 
written  in  the  following  form 

fie  =  0  -  tan^UJU,).  (2.1) 

In  this  study,  a  cross  section  on  the  stator  at  55%  radius  of  rotor  has  been  selected 
for  both  surface  pressure  and  velocity  measurements  (Fig.  2-2(a)).  The  zero-lift-line 
angle,  0,  at  rotor's  55%R  section  is  approximated  as  26° .  Meanwhile,  the  corresponding 
chord  length  of  rotor  at  55 %R  is  measured  to  be  1.8  cm. 
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Gap  between  Rotor  and  Stator 

The  gap  between  the  rotor  and  stator  is  fixed  in  all  investigated  cases  but  changes 
along  the  rotor’s  spanwise  direction  due  to  the  variation  in  both  chord-length  and 
incidence  of  the  rotor’s  blades.  Approximately,  the  gap  varied  from  0.5  cm  to  0.865  cm 
along  the  direction  from  the  rotor's  root  to  its  tip.  Hence,  a  constant  distance,  instead  of 
the  gap,  is  measured  to  be  1  cm  between  the  center-chord-line  along  the  rotor’s  span  and 
the  leading  edge  of  stator  blade,  as  shown  in  Fig.2-3.  The  gap  at  55%R  of  rotor  is 
measured  as  0.71  cm.  Thus,  this  gap  is  39%  of  the  corresponding  rotor’s  chord  length 
(1.8  cm)  and  14%  of  the  stator’s  chord  length  (5cm). 

It  is  noticed  that,  from  the  setup  sketched  in  Fig.2-3,  this  rotor  and  stator 
arrangement  is  a  typical  stage  in  an  axial  compressor  (Turton  1984). 


2.2.  Measurements 

Three  types  of  measurements  have  been  conducted  in  this  research.  They  are 
downstream  noise  measurement,  unsteady  surface  pressure  measurement  on  the  stator, 
and  velocity  measurement  around  the  stator. 
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2.2.1.  Downstream  Noise  Measurement 

In  order  to  identify  the  downstream  noise  signature,  the  microphone 
measurements  are  set  up  at  two  different  pipe-sections  downstream  from  the  rotor-stator 
stage,  as  illustrated  in  Fig.  2-2(c).  First  downstream  measurement  is  accomplished  by  a 
quarter-inch-diameter  Bruel  &  Kjaer  (B&K)  type-4135  condenser  microphone  flush- 
mounted  on  the  pipe  wall  which  was  0.7  m  distant  from  the  stator’s  leading  edge.  In 
addition,  two  quarter-inch-diameter  B&K  (Type  4135  and  Type  4136)  condenser 
microphones  are  placed  90°  apart  and  also  flush-mounted  on  the  pipe  wall  which  is  0.23 
m  downstream  from  the  leading  edge  of  stator. 

2.2.2.  Surface  Pressure  Measurement 

A  total  of  37  pressure  holes  on  both  sides  of  the  stator  blades  are  used  for  the 
surface  pressure  measurement:  19  pressure  holes  on  the  pressure  sides  and  18  on  the 
suction  sides.  They  are  located  at  the  mid-span  of  a  single  stator  blade  which  was  55% 
radius  of  rotor  distant  from  the  center  line  of  the  test  section.  These  pressure  holes  are 
drilled  to  meet  channels  underneath  the  stator  blade,  and  these  channels  terminate  up 
with  connections  to  metal  tubes  at  the  end  side  of  the  stator  blade  (Fig.  2-4(a)).  The 
pressure  holes  are  typically  in  two  different  diameters:  0.254  mm  (0.01  inch)  around  the 
leading  edge  and  0.508  mm  (0.02  inch)  otherwise  on  both  sides.  Due  to  the  difficulty  of 
drilling  the  channels  near  the  trailing  edge,  it  was  unable  to  measure  the  surface  pressure 
at  the  rear  15%C  on  both  sides  of  stator  (Fig.  2-4(a)).  Fortunately,  this  deficiency  does 
not  influence  the  current  study  since  the  flow  interaction  happens  mainly  around  the 
stator’s  leading-edge  portion,  as  we  will  see  in  later  chapters. 

By  connecting  to  each  pressure  hole,  the  surface  pressure  measurement  on  the 
stator  consists  of  two  separate  parts:  1)  mean  pressure  measured  by  a  static-pressure  tube 
with  a  1-torr-range,  differential-type  MKS  pressure  transducer,  and  2)  pressure 
fluctuation  measured  by  a  quarter-inch-diameter  B&K  type-4136  condenser  microphone 
via  a  transmission  line.  The  transmission  line  is  composed  of  the  surface  channel  in  the 
stator  blade  from  the  pressure  hole  to  the  end-side  metal  tube  and  a  Tygon  tube 
(1/16"I.D.,  9.7  cm  in  length)  connecting  the  metal  tube  to  a  well-sealed  microphone 
adapter  for  the  B&K  condenser  microphone,  in  Fig.  2-4(a). 

Before  taking  surface  pressure  measurements,  a  calibration  process  for  the 
transmission  line  (Fig.2-4(a))  should  be  conducted  in  advance.  This  is  due  to  energy  loss 
and  phase  lag  of  the  pressure  wave  propagating  through  the  transmission  line.  The 
pressure-fluctuation  measurements  by  the  condenser  microphone  should  be  compensated 
in  order  to  recover  the  distorted  signal  back  to  the  actual  signal  on  the  stator  surface. 
Accordingly,  gain  and  phase  characteristics  or  frequency  responses  for  the  transmission 
line  should  be  determined  prior  to  the  compensating  process.  There  are  37  pressure  holes 
used  for  surface  pressure  measurement;  hence,  37  sets  of  gain  and  phase  responses 
should  be  compensated  in  order  to  yield  the  correct  pressure  distributions  on  the  stator 
surface. 
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Experimental  Procedure  for  Pressure- Fluctuation  Compensation 

A  conceptual  description  to  the  calibration  process  of  obtaining  the  frequency 
responses  for  a  transmission  line  is  illustrated  in  Fig.  2-4(b).  In  the  real  setup,  a  7"D  sub¬ 
woofer  is  chosen  for  the  sound  source  and  attaches  to  one  end  of  the  7"I.D.  plexiglass 
pipe.  It  generates  sinusoidal  waves  in  a  frequency  range  from  35  to  3000  Hz  and  is 
driven  by  a  wave  generator.  The  blade  with  a  transmission  line,  as  illustrated  in  Fig.2- 
4(a),  and  a  reference  microphone  are  both  placed  at  the  same  pipe  section  which  is  19  cm 
downstream  from  the  sub-woofer.  On  the  plane  of  this  pipe  section,  a  plane  acoustic 
wave  has  been  detected  and  checked  by  the  reference  microphone.  Both  microphone 
signals  were  recorded  via  analog-to-digital  converter  to  the  computer.  By  varying  the 
frequencies,  we  are  able  to  measure  both  reference  and  testing  signals  from  two 
microphones  for  determination  of  frequency  responses  of  this  transmission  line.  The 
procedure  to  obtain  the  frequency  responses  will  be  described  in  Section  2.3.2.  A  single¬ 
input/single-output  model  is  illustrated  to  describe  the  transmission-line  system  in  Fig.  2- 
4(c).  Fig.  2-4(d)  schematically  illustrates  a  linear-system-identification  process  by  a 
finite-duration-impulse-response  (FIR)  filter  for  compensating  the  distorted  pressure- 
fluctuation  signals,  which  will  be  discussed  in  detail  later  in  Section  2.3.3. 
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2.2.3.  Velocity  Measurement 

A  special,  single-component  hot-wire  probe,  with  broaches  bent  45°  from  their 
original,  straight  positions,  has  been  utilized  in  the  velocity  measurements  around  the 
stator,  as  shown  in  Fig.2-2(a)  &  2-5.  This  special  probe  expedited  velocity  measurement, 
especially  around  the  region  adjacent  to  the  leading  edge  of  stator.  In  this  study,  the  wire 
between  two  broaches  on  the  probe  was  set  parallel  to  the  spanwise  direction  of  the  stator 
(Fig.  2-5).  When  the  hot-wire  probe  was  moved  close  to  the  stator  surface,  this  setup 
provided  a  better  evaluation  of  the  streamwise  velocity  component  because  of  the 
physical  weakness  of  the  normal  velocity  component  close  to  a  solid  boundary  and  the 
poor  resolution  of  the  spanwise  velocity  component  due  to  the  worst  flow  heat-transfer 
effect  along  the  wire  (Blackwelder  1981). 

A  three-axis  traverse  system  was  used  to  support  and  move  the  hot-wire  probe, 
and  it  was  driven  by  PC  through  a  three-axis  controller  manufactured  by  Anaheim 
Automation,  Inc. 

Since  the  stator  blades  possessed  both  a  constant  NACAOO  16-cross-section  and 
constant  incidence  along  its  span,  we  may  select  cross  sections  away  from  both  the  hub 
and  the  tunnel  wall  to  investigate  the  unsteady  flow  field  mainly  induced  by  rotor-stator 
interaction.  The  velocity  measurements  were  taken  specifically  around  the  front  30% 
chord  length  of  stator’s  blade  at  55%  radius  of  rotor.  The  measurements  were 
accomplished  by  generating  a  fine  mesh  of  approximately  4500  hot-wire  stations  which 
will  be  shown  in  the  later  chapters. 

2.3.  Data  Processing  Technique 

The  data  acquisition  system  includes  a  PC  486DX50  and  a  Dell  Pentium  90  with 
a  12-bit,  16-channel  analog-to-digital  board  manufactured  by  RC  Electronics,  Inc.  The 
maximum  sampling  frequency  is  1  MHz.  This  maximum  sampling  frequency  is  sufficient 
enough  for  sampling  data  in  this  study  because  the  maximum  value  of  blade  passing 
frequency  (BPF)  will  not  exceed  133Hz  (8000  RPM)  for  rotor  RPM  less  than  4000. 

2.3.1.  Phase-Averaging  Technique 

As  mentioned  previously,  a  square- wave  reference  related  to  the  rotor’s  motion 
can  be  employed  for  data  analysis  with  a  phase- averaging  technique.  Typical  square- 
wave  reference  and  measurement  signals  in  this  study  are  shown  in  Fig.2-6(a). 

Phase-Averaging  Algorithm 

The  phase-averaging  technique  is  introduced  for  extracting  the  deterministic  part 
of  the  signal.  Basically,  it  has  a  resemblance  to  the  ensemble  averaging  method  except 
that  the  phase  average  is  based  on  an  intrinsic  period,  for  example,  T  in  Fig.2-6(a). 

The  phase-averaging  algorithm  is  defined  as,  (Fujita  and  Kovasznay  1974), 
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where  is  the  period,  and  N  is  the  number  of  cycles  or  samples  for  phase-averaging.  In 
practice,  the  choice  of  cycle  number  depends  on  a  compromise  between  the  accuracy  of 
the  result  and  the  time  length  available  to  accomplish  one  experiment.  Generally 
speaking,  500  <N  <  1000  is  the  typical  range  of  cycle  number  chosen.  Fig.  2-6(b)  shows 
the  phase-averaging  results  of  the  signals  shown  in  Fig.  2-6(a),  by  using  approximately 
600  cycles  to  smooth  out  the  random  fluctuation.  At  first  sight  of  the  raw  signals,  it  is 
noticed  that  the  noise  shown  in  part  (a)  of  Fig.  2-6  has  no  significantly  distinguishable 
repeating  pattern  as  the  other  signals.  However,  in  Fig.  2-6(b),  a  high  peak  in  the  phase- 
averaged  noise  signature  emerges  from  the  random  fluctuations  which  are  mostly  the 
pressure  fluctuations  due  to  small  eddies  passing  by  the  flush-mounted  microphones 
(Lighthill  1962,  Goldstein  1983).  The  random  fluctuations  are  filtered  by  phase¬ 
averaging  technique. 

In  the  present  study,  150,  500  and  1000  cycles  were  adopted  for  phase-averaging 
the  velocity,  surface  pressure  and  noise  signals,  respectively. 


2.3.2.  Frequency  Responses  of  Transmission  Line 

In  Section  2.2.2,  we  learn  that,  in  order  to  compensate  the  distorted  pressure- 
fluctuation  signals,  the  frequency  responses  of  transmission  lines  (Fig.  2-4(a))  should  be 
determined  in  advance.  Both  the  reference-  and  testing-microphone  signals  (Fig.  2-4(b) 
to  (d))  were  used  to  obtain  both  the  gain  and  phase  frequency  responses. 

The  gain  response  at  a  certain  frequency  can  be  obtained  by  taking  the  ratio  of 
both  microphone  signal  magnitudes  which  are  calculated  from  their  energy  spectra  by  a 
fast  Fourier  transform  (FFT)  technique  (Brigham  1988,  Proakis  and  Manolakis  1992). 
Meanwhile,  the  phase  lag  or  response  at  the  same  testing  frequency  is  acquired  by 
directly  computing  the  phase  angle  of  the  one-sided  cross-spectral  density  function  from 
both  signals  (Bendat  and  Piersol  1986). 

The  gain  and  phase  responses  for  two  typical  pressure  holes  (diam. =0.254  mm 
and  0.508  mm)  are  plotted  in  Fig.  2-7. 


2.3.3.  Compensating  Procedure  for  Pressure  Fluctuation 

Three  processes  are  needed  to  accomplish  the  compensation  of  transmission 
distortion:  (  I  )  the  linear-inverse-system  identification  process  for  a  transmission  line 
(Proakis  and  Manolakis  1992),  which  uses  a  digital  finite-duration  impulse  response 
(FIR)  filter  to  identify  a  linear  system  as  an  inverse  system  of  the  gain  response  of  a 
transmission  line,  ( II )  the  gain- compensating  process  of  pressure  fluctuation  due  to  both 
energy  loss  and  phase  lag  in  the  transmission  line,  which  mainly  applies  a  convolution 
process  between  the  raw  signal  of  pressure  fluctuation  and  the  compensating  FIR  filter 
obtained  in  PROCESS(  I  ),  and  (  III  )  the  phase-compensating  process  of  pressure 
fluctuation  due  to  phase  lag  in  the  transmission  line,  which  transforms  the  gain- 
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compensated  signal  of  the  pressure  fluctuation  to  the  frequency  domain  by  Fourier 
transform,  after  PROCESS(  II ),  in  order  to  subtract  the  phase  lag  at  each  frequency,  and 
then  obtain  the  well-compensated  signal  of  pressure  fluctuation  in  the  time  domain  by 
inverse  Fourier  transform. 

With  knowledge  of  the  frequency  responses  for  all  transmission  lines,  we  may 
start  to  drive  the  compensating  filters  for  each  transmission  line.  A  digital  finite-duration 
impulse  response,  FIR,  filter  (McClellan  et  al.  1975,  Rabiner  et  al.  1975,  Rabiner  and 
Gold  1975)  is  employed  as  a  linear,  time-invariant  (LTI)  system  to  accomplish  the 
linear-inverse-system  identification  process  for  a  transmission  line  (Fig.  2-4(d)).  The 
major  advantage  of  adopting  the  digital  FIR  filter  is  that  it  possesses  exact  linearity  of 
phase.  This  makes  both  the  gain-  and  phase-compensating  processes  easier  and  more 
feasible,  specifically  in  the  phase  compensation  part.  It  should  be  noticed  that  a  low-pass 
filter  used  in  this  procedure,  as  will  be  described  later  in  this  section,  was  also  designed 
by  the  FIR  filter  technique. 

The  convolution  technique  is  used  for  processing  the  raw  pressure-fluctuation 
signal  via  FIR  filters,  including  the  low-pass  FIR  filter,  described  in  the  preceding  sub¬ 
section.  Next,  by  transforming  the  gain-compensated  signal  with  a  fast-Fourier- 
transform  (FFT)  technique,  the  fully  compensated  signal  was  recovered  by  directly 
subtracting  the  phase  lag  at  each  frequency  with  the  phase  response  information,  and 
could  be  obtained  by  using  an  inverse-fast- Fourier- transform  (IFFT)  technique  (Brigham 
1988,  Proakis  and  Manolakis,  1992). 

Steps  of  the  Compensating  Procedure 

The  steps  of  the  compensating  procedure  for  pressure-fluctuation  signals  are 
listed  as  follows: 

(1)  find  the  digital  compensating  filter  for  each  transmission  line  within  the 
desired  frequency  range, 

(2)  start  the  gain  compensating  for  each  pressure-fluctuation  signal  obtained 
from  each  transmission  line  by  taking  the  convolution  of  the  signal  and  its 
corresponding  digital  compensating  filter, 

(3)  use  a  low-pass  FIR  filter  to  keep  the  gain-compensated  signals  within  the 
desired  frequency  range, 

(4)  multiply  an  appropriate  constant  to  obtain  the  correct  gain  for  signals 
obtained  in  STEP  (3), 

(5)  take  the  FFT  of  the  gain-compensated  single  to  obtain  both  amplitude  and 
phase  angle  spectra  in  the  frequency  domain, 

(6)  get  the  phase  angle  at  each  frequency  and  subtract  the  corresponding 
phase  lag  measured  by  the  same  transmission  line,  and 

(7)  take  IFFT  of  both  amplitude  and  phase-compensated  spectra  to  yield  the 
fully-compensated  pressure-fluctuation  signals  in  time  domain. 

In  this  procedure,  STEP(l)  is  the  linear-inverse-system  identification  process, 
using  a  FIR  filter,  for  a  transmission  line.  STEP(2)  to  (4)  are  related  to  the  gain- 
compensating  process  of  pressure  fluctuation  due  to  the  frequency  loss  in  transmission 
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lines,  while  STEP(5)  to  (7)  are  the  phase-compensating  process  for  pressure  fluctuation. 
The  gain  responses  in  Fig.  2-7  (a),  for  two  typical  pressure  holes  in  this  study,  are 
compensated  right  after  STEP(4)  and  plotted  in  Fig.  2-8.  The  flat  lines  in  Fig.  2-8(a)  & 
(b)  represent  the  well-compensated  gain  responses  which  were  obtained  by  multiplying 
the  gain  responses  in  Fig.  2-7 (a)  with  the  linear-inverse-system  function  by  FIR  filters. 
An  example  used  to  demonstrate  the  compensating  procedure  can  be  found  in  Hu  and 
Ho’s  report  (1997). 
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CHAPTER  3 
NOISE  MEASUREMENT 


In  Section  2.2.1,  we  have  described  that  the  microphone  measurements  were 
conducted  for  identifying  the  sound  wave  at  two  pipe  sections  which  were  0.23m  and 
0.7m  downstream  from  the  stator’s  leading  edge.  Since  the  plane  pressure  wave  measured 
by  microphones  at  the  0.23m  pipe  section  did  not  propagate  at  the  acoustic  speed  (Hu  and 
Ho  1997),  hence  it  can  not  be  recognized  as  the  sound  wave.  Meanwhile,  due  to 
identifying  a  plane  acoustic  wave  arriving  at  the  0.7m  pipe  section  (Hu  and  Ho  1997), 
only  the  noise  signatures  were  measured  by  Microphone  A  at  position  #  3  of  this  pipe 
section  (Fig.  2-2(b)). 

Delay  Time,  R/a^,  and  phase  advancing 

Since  these  noise  signatures  are  measured  at  the  pipe  section  0.7m  downstream 
from  the  stator’s  leading  edge,  the  delay  times,  R/a^,  for  this  pipe  section  (as  indicated  in 
Chapter  1)  is  calculated  to  be  2.06  msec,  based  on  the  value  of  sound  speed  of  340  m/s  in 
air  at  20°C.  This  delay  time  are  needed  to  advance  the  measured  noise  signals  in  order  to 
be  correlated  with  the  unsteady  flow  field  around  the  rotor-stator  stage.  Notice  that  the 
noise  signatures  in  this  study  will  be  plotted  against  the  dimensionless  axis  of  phase  angle 
</>.  Thus,  we  have  to  advance  this  delay  time  in  phase  angle  ^ ,  i.e.  a  normalized  time.  The 
phase  advancing  is  calculated  by  {delay  time!  T  )-360°,  where  is  defined  as  the  period 
based  on  BPF.  The  noise  signatures  after  the  phase  advancing  are  shown  in  Fig.3-1 
through  3-4. 

In  this  chapter,  we  will  define  an  active  phase  interval  by  comparing  noise 
signatures  measured  with  and  without  the  presence  of  the  stator.  This  active  phase  interval 
will  be  utilized  to  correlate  the  acoustic  field  to  the  unsteady  flow  field  in  the  later 
chapters. 

3.1.  Operating  Conditions  for  Noise  Measurements 

In  this  investigation,  there  are  four  freestream  speeds,  five  rotor  speeds  at  55%R, 
and  six  different  stator’s  AoA  chosen  in  the  downstream  noise  measurement,  and  therefore 


Table  3-1. 


@  55%R 
of  rotor  (m/s) 


Pe  and  Ut  at  55%R  of  rotor  for  noise  measurements. 


Uo 

(m/s) 

4.66 

5.62 

8.43 

11.24 

{],  @  55%R 
of  rotor  (m/s) 


U,  (m/s) 


6.11 

-11.33° 

7.68 

-16.61° 

8.30 

-28.07° 

10.41 

-35.47° 

12.79 

7.65 

-  5.35° 

8.96 

-10.30° 

9.49 

-21.78° 

11.38 

-29.76° 

13.60 

12.23 

5.14° 

13.09 

1.32° 

13.46 

-  8.58° 

14.85 

-16.58° 

16.61 

15.29 

9.05° 

15.98 

5.82° 

16.29 

-  2.87° 

17.46 

-10.32° 

18.98 

19.11 

12.30° 

19.67 

9.61° 

19.92 

2.20° 

20.89 

-  4.46° 

22.17 

we  have  one-hundred-and-twenty  operating  conditions  set  up  for  studying  the  interaction- 
noise  signatures  in  a  range  of  -35.5°  <Pe  <  12.5°  and  -15°  <  Oi  <  10°. 

The  freestream  speeds  and  the  rotor  speeds  are  shown  in  Table  3-1  8l  Fig.  3-5. 
The  six  stator’s  AoA  as  are  0°,  +5°,  ±10°  and  -15°,  respectively. 

3.2.  Noise  Signatures 

After  being  phase-averaged  and  advanced  with  the  delay  time,  the  noise  signatures 
are  shown  in  Figures  3-1,  3-2,  3-3  and  3-4.  In  these  figures,  the  horizontal  axis  is  a  non- 
dimensionalized  time  axis  in  phase  angle,  (f),  while  the  vertical  axis  represents  noise 
magnitude  with  a  dimensional  unit.  Pa.  The  phase  angle,  (j) ,  shown  in  Fig.4-3,  is 
normalized  by  the  individual  T  value  of  each  case  and  re-scaled  from  -180°  to  180°.  The 
phase  angle  ^  at  0°  is  defined  as  when  the  rotor’s  trailing-edge  position  is  aligned  with 
the  stator  leading  edge  which  is  at  0°  AoA  (Fig.2-2(b)  and  Fig.2-7(b)).  Similarly,  (f>  at 
±180°  is  defined  as  rotor's  position  perpendicular  to  the  stator,  as  shown  in  Fig.  2-7(b). 

3.2.1.  Noise  Peak 

By  carefully  examining  the  noise  signatures  in  Figures  3-1  through  3-4,  a 
prominent  noise  peak  is  consistently  observed  before  the  zero-degree  phase  angle  for 
cases  where  Pe  >  -22°  with  all  six  stator's  AoA.  This  noise  peak  is  an  indicator  for  both 
event  time  and  strength  of  noise  generation  due  to  rotor-stator  interaction. 

Alignment  of  Noise  Peaks  for  Cases:  p^  >  -22  ° 

The  first  observation,  from  noise  signatures  in  Fig.  3- 1  to  3-4,  is  an  alignment  of 
noise  peaks  at  about  -30°  phase  angle  in  all  cases  where  p^  >  -22°.  In  other  words,  the 
time  when  the  peak  noise  happens  is  irrelevant  to  Vo ,  Vr  and  P^,  for  all  cases  where  P^  >  - 
22°.  This  result  is  probably  due  to  both  the  similar  flow  fields  and  the  fixed  gap  between 
the  rotor  and  the  stator  in  this  study.  For  those  cases  where  p^  <  -22°,  the  unsteady  flow 
fields  corresponding  to  the  interaction-noise  generation  might  be  fairly  different  from 
those  in  p^  >  -22°.  Accordingly,  in  this  study,  the  operating  conditions  for  -22°  <  P^<  12° 
is  used  for  a  comprehensive  survey  of  interaction-noise  signatures. 


16 


Noise  Peaks  vs.  Velocities  for  Cases:  /S,  >  -22  ° 

Our  next  interest  is  to  look  for  a  relationship  between  the  noise-peak  value  and 
the  velocities.  In  general,  for  the  cases  where  fie  >  -22°,  the  peak-noise  value  goes  up  as 
long  as  one  or  both  of  Ug  and  Ur  increase  under  a  fixed  stator's  AoA,  as,  by  looking  into 
the  noise  signatures  in  Fig.  3-1  to  3-4.  Although  the  event  time  of  the  noise  peaks  appear 
irrelevant  to  both  Uo  and  Ur  in  these  cases,  we  still  expect  a  relationship  between  noise- 
peak  values  and  the  velocity  magnitudes.  A  simple  relation  between  peak  values  and 
velocities  can  be  found  by  a  power-fitting  method  as  follows; 

Y  =  aX^,  or  \og(Y)  =  b  log(  +  A  ,  where  A  =  log(a).  (3.1) 

Choosing  7  for  the  peak  value,  the  coefficients  a’s  and  power  i’s  are  listed  in  Table  3-2 
with  X  for  freestream  velocity  Uo,  rotor’s  velocity  Ur  at  55%R,  and  the  resultant  velocity 
Ut  at  55%R,  respectively,  under  different  stator’s  AoA,  a,. 


X 

II 

X  =  U,@  55%R 

X  =  U.@  55%R 

a,  * 

a 

b 

a 

b 

a 

b 

-15° 

0.0163 

2.031 

0.0098 

1.739 

0.0015 

2.339 

-10° 

0.0210 

1.828 

0.0077 

1.787 

0.0013 

2.347 

-5° 

0.0146 

1.941 

0.0090 

1.664 

0.0014 

2.263 

0° 

0.0141 

1.975 

0.0414 

1.063 

0.0077 

1.629 

5° 

0.0077 

2.229 

0.0234 

1.241 

0.0038 

1.853 

10° 

0.0072 

2.195 

0.0208 

1.237 

0.0035 

1.829 

*  a, :  the  stator’s  AoA. 


From  Table  3-2,  it  shows  that  the  power  Ns  have  average  values  of  2.03,  1.43  and 
2.03  for  velocities  Uo,  Ur  at  55%R,  Ut  at  55%R,  respectively.  This  means  that  noise-peak 
value  increases  as  any  of  the  velocities  increases.  Since  (/,  is  the  resultant  velocity  of  both 
Uo  and  Ur,  hence  we  may  assert  that  it  directly  contributes  to  the  noise-peak  magnitude. 
Fig.  3-6  shows  the  relationship  between  peak  value  and  C/,.  Thus,  the  peak  value  is 
proportional  to  t//,  i.e.. 

Peak  Value  oc  f//  (3.2a) 


or. 


Peak  Value  oc  q, , 


(3.2b) 


where  qt  is  a  local  dynamic  pressure  defined  as  .  Physically,  it  follows  that  the 

peak  noise  is  directly  related  to  the  local  flow  field,  i.e.  qt,  due  to  rotor-stator  interaction. 


Table  3-3.  Y  =  (Peak  Value /q,)  &X-  Velocities  in  Eg. (3.1)  in  cases;  fie  >  -22 


X  =  U„ 


U  @  55%R 


X  =  U.@  55%R 
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a,  * 

a 

b 

a 

b 

a 

b 

-15° 

7.66x10^ 

1.087 

4.47x10"' 

0.116 

2.23x10"' 

0.370 

-10° 

1.16x10' 

0.817 

3.66x10' 

0.156 

2.10x10"' 

0.357 

-5° 

6.62x10^ 

1.012 

4.35x10"' 

0.011 

2.22x10"' 

0.265 

0° 

5.63x10^ 

1.110 

1.90x10"' 

-0.570 

1.13x10"' 

-0.338 

5° 

5.11x10^ 

1.128 

1.10x10"' 

-0.376 

6.34x10"' 

-0.146 

10° 

3.81x10-' 

1.196 

9.17x10"' 

-0.368 

5.31x10"' 

-0.140 

*  a, :  the  stator’s  AoA. 


If  we  normalize  the  peak  values  with  their  corresponding  ^,’s  for  the  cases  where 
>  -22°,  then  we  can  apply  the  power-fitting  equation  (3.1)  again  to  find  the 
relationship  between  normalized  peak  values  and  velocities.  The  corresponding 
coefficient  a’s  and  power  b's  are  listed  in  Table  3-3. 

From  Table  3-3,  average  values  of  power  i’s  for  Uo ,  Ur  and  C/,  are  1.06,  -0.17  and 
0.06,  respectively.  It  is  apparently  that  the  freestream  velocity  Uo  remains  linearly 
proportional  to  the  noise-peak  values  at  different  stator’s  AoA.  This  power  relationship 
between  normalized  peak  value  and  freestream  velocity  Uo  at  different  stator’s  AoA  is 
plotted  in  Fig. 3 -7,  and  can  be  expressed  as 

(Peak  Value  /q,)cc  Uo  (3.3  a) 

or 

Peak  Value  oc  q,  Uo  (3.3b) 

Before  discussing  the  physical  significance  of  Eq.  (3.3),  we  should  examine  if  there  still 
exists  any  power  dependency  between  peak  values  and  velocities.  Once  again,  the 
coefficient  a’s  and  power  b’s  for  Y  =  (Peak  Value /q,  Uo)  are  listed  in  Table  3-4. 


Table  3-4.  Y  =  (Peak  Value/ qr  Up)  &  X  =  Velocities  in  Eq.(3.1)  in  cases;  fie  >  -22°. 


X 

II 

X  =  U 

@  55%R 

X  =  U 

@  55%R 

a,  * 

a 

b 

a 

b 

a 

b 

-15° 

0.0010 

-3.06x10"" 

0.0010 

-6.54x10" 

0.0010 

-1.49x10"" 

-10° 

0.0010 

-3.06x10" 

0.0010 

-6.54x10"" 

0.0010 

-1.49x10" 

-5° 

0.0010 

-3.06x10" 

0.0010 

-6.54x10" 

0.0010 

-1.49x10" 

0° 

0.0014 

-0.142 

0.0014 

-0.121 

0.0017 

-0.179 

5° 

0.0010 

1.79x10"" 

0.0010 

-2.39x10" 

0.0010 

-1.27x10" 

10° 

0.0010 

-2.15x10" 

0.0010 

-1.69x10" 

0.0010 

-4.58x10" 

*  a, ;  the  stator’s  AoA. 


From  Table  3-4,  we  observe  that  all  the  power  b’s  approach  zero  for  three 
velocities  except  for  =  0°.  In  other  words,  Eq.(3.3b)  shows  a  general  power  relationship 
between  noise-peak  magnitude  and  velocities  for  cases  where  >  -22°.  We  may  rewrite 
Eq.(3.3)  as 
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Peak  Value  -  const,  q,  Uo  (3.4) 

In  general,  Eqs.(3.3)  and  (3.4)  not  only  describe  an  amplitude  relationship  between 
noise-peak  values  and  velocities  but  also  provide  a  pressure  parameter,  q,,  for  non- 
dimensionalizing  the  noise  signatures  for  those  cases  where  >  -22°  in  Fig.  3-1,  3-2,  3-3 
and  3-4.  In  addition,  from  both  equations,  it  also  shows  that  these  dimensionless  peak 
amplitudes  are  still  proportional  to  the  freestream  velocity  Uo.  In  other  words,  this  implies 
that  the  dimension  of  the  proportional  constant  in  Eq.(3.4)  has  a  unit  as  an  inverse  of  a 
velocity  scale.  Furthermore,  from  Eq.(3.4),  we  observe  that  the  contribution  to  the  noise- 
peak  values  from  the  rotor’s  velocity  Ur  is  implicitly  included  in  qt  which  is  equal  to 
\p,U^  or\p,{Ul+U^,). 


Noise  Peak  vs.  for  Cases:  fie  >  -22  ° 

Since  fie  is  directly  related  to  the  velocity  ratio  Uo/Ur  and  defined,  in  Eq.(2. 1),  as 

Pe  =  e-tan-\UJU^, 


or 

UolUr  =  ton(e-pj,  (3.5) 

we  can  derive,  from  Eq.(3.3),  a  form  that  relates  the  peak  value  to  Pe  in  the  following. 


For  Ut  is  the  resultant  velocity  of  both  Uo  and  Ur,  then  we  can  rewrite  Eq.(3.3b)  as 

Peak  Value  <x:  q,Uo  =  2  Po  Uo 

=  \po{Uo'+Ur^)Uo 
=  kPoUo\\+Ur'lUo')Uo. 

By  substituting  Eq.(3.5)  for  the  velocity  ratio  UrIUo ,  we  obtain 

Peak  Value  ex  (70(1+ - — — —)Uo  (3.6) 

tan  (0-pJ 


where  qo-  2  PoUo^  is  the  freestream  dynamic  pressure.  Comparing  Eq.(3.6)  and  Eq.(3.3), 
we  have 


% 


=  qo{\+ 


1 

to«^(0-P,) 


). 


(3.7) 


From  Eq.(3.6)  and  (3.7),  it  turns  out  that  the  contribution  to  the  noise-peak  value 
from  the  rotor’s  effective  AoA  Pe  is  implicitly  included  in  the  local  dynamic  pressure,  qt, 
like  the  rotor’s  velocity,  Ur.  Since,  in  the  present  study,  Pe  at  55%R  ranges  from  -22°  to 
12°,  with  a  fixed  0  of  26°  at  the  same  radius,  the  value  in  the  parentheses  on  the  right- 
hand  side  of  Eq.(3.6)  varies  from  1.8  to  17.1,  respectively.  Apparently,  from  Eq.(3.6),  the 
noise-peak  magnitude  increases  as  Pe  increases  with  a  constant  6 ,  for  all  cases  where  Pe 
>  -22°. 
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Peak  Noise  vs.  for  Cases:  0°  <  p^<  10° 

One  interesting  thing  observed  from  noise  signatures  in  Fig.  3-1  to  3-4  is  that,  for 
cases  that  0°  <  /?^  <  10°,  the  noise  peak  increases  as  the  stator’s  AoA  decreases. 
However,  there  does  not  exist  a  similar  relationship  between  the  peak  values  and  a  in 
these  cases  as  Eq.(3.6).  Nevertheless,  this  observed  phenomenon  still  reveals  that  a  is  an 
important  parameter  to  the  noise  generation  for  cases  where  0°  <  <  10°.  Further 

investigation  of  the  unsteady  flow  field  induced  by  rotor-stator  interaction  is  necessary  to 
understand  this  phenomenon. 

3.2.2.  Interaction  Period  for  Cases:  Pe  >  -22°. 

As  indicated  previously,  the  noise  peak  represents  a  maximum  noise  generated  by 
rotor-stator  interaction  for  cases  where  Pe  >  -22°.  Its  existence  indicates  a  reference 
phase  angle  for  the  active  interaction  period.  Since  most  noise  peaks,  for  the  cases  where 
Pe  >  -22°,  lie  in  an  interval  of  phase  angle  -27°±3°,  we  may,  for  the  sake  of  simplicity, 
choose  a  phase  angle  -30°  as  the  upper  bound  of  the  active  interaction-noise  period. 
Meanwhile,  a  sudden  rise  exists  commonly  at  ^  =  -60°.  We  choose  this  phase  angle  as  the 
lower  bound  of  the  interaction  period.  We,  hence,  define  an  interaction  period  from  a 
phase  angle  -60°  to  -30°  for  the  remainder  of  this  study. 

3.2.3,  Two  Typical  Noise  Patterns  in  Cases:  Pe  >  -22°. 

In  this  subsection,  we  are  going  to  compare  two  typical  noise  patterns  observed 
from  cases  where  pe  >  -22°.  Two  representative  noise  signatures  are  plotted  in  Fig.  3-8. 

These  two  signatures  are  normalized  by  their  corresponding  local  dynamic 
pressure  qt  and  chosen  from  cases  where  Pe  =  -10.30°  in  Fig.  3-2(b)  and  Pe  =  9.61°  in 
Fig.  3-2(b)  with  a  fixed  as  =  10°.  Simply,  we  use  Pe  =  ±10°  to  distinguish  these  two  cases 
for  discussion. 

In  Fig.  3-8,  we  observe  that  these  two  normalized  noise  signatures  are  apparently 
different  in  both  magnitude  and  pattern.  The  dimensonless  noise  magnitude  in  the  case  Pe 
-  -10°  is  obviously  higher  than  the  one  in  the  case  Pe  =  10°.  For  the  case  Pe  =  -10°,  the 
noise  starts  decreasing  before  phase  angle  ^  =  -90°,  reaches  a  valley  at  ^  =  -60°,  and  then 
rises  abruptly  to  a  peak  zX  <f>  =  -30°.  However,  for  the  case  Pe  =  10°,  the  noise  rises 
gradually  from  a  valley  at  ^  =  -90°,  and  then,  after  reaching  a  higher-noise  level  at  ^  =  - 
48°,  increases  rapidly  to  a  peak  at  ^  =  -30°.  Differences  observed  in  both  the  magnitude 
and  the  pattern  suggest  that  there  may  be  two  distinct  mechanisms  which  generate  the 
interaction  noise  during  the  rotor’s  passage  in  the  vicinity  of  the  stator’s  leading  edge. 
These  two  typical  noise  signatures  generally  exist  in  two  categories  where  Pe  <  0°  and  Pe 
>  0°,  respectively.  Further  discussion  will  be  postponed  until  investigation  of  the 
interaction-flow  fields  in  the  later  chapters. 
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3.3.  Active  Interaction  Phase  Interval 


In  sub-section  3.2.2,  we  have  already  defined  the  interaction  period,  for  Pe>  - 
22°,  by  locating  a  sudden  rise  of  peak  noise  from  -60°  to  -30°  in  phase  angle  (j).  In  order 
to  confirm  the  relationship  of  this  peak  noise  to  the  rotor-stator  interaction  noise,  we  have 
to  examine  noise  signatures  which  are  measured  with  and  without  the  presence  of  a  stator 
for  two  cases:  /3e  =  -10.30°  and  5.82°. 

Comparison  of  noise  signatures  with  and  without  Stator 

In  Fig.3-9(a)  and  (b),  it  is  obvious  that,  without  the  stator,  the  peak  noise  does  not 
exist  for  both  cases.  Accordingly,  from  Fig.3-9,  the  peak  noise,  with  the  presence  of 
stator,  is  identified  as  the  rotor-stator  interaction  noise. 

Active  interaction-noise  phase  interval  for  Pe  >  -22  ° 

As  a  consequence,  from  Fig.3-9,  the  interaction  noise  arises  abruptly  from  a 
valley  to  a  peak  before  the  zero-degree  phase  angle,  when  compared  to  the  noise 
measured  without  a  stator.  Therefore,  the  active  time  interval  of  noise  interaction  or  the 
active  interaction  phase  interval  is  defined  between  the  corresponding  phase  angles  of 
these  two  values,  and  can  be  found  in  Fig.  3-9.  Therefore,  from  Fig.  3*9(a)  &  (b),  the 
active  interaction  phase  interval  for  both  cases  ranges  from  -60°  to  -30°  in  phase  angle  (}>. 
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3.4.  Summary 

According  to  the  discussion  in  previous  sections,  we  may  summarize  as  follows. 
First,  it  is  noted  that  the  active  interaction  phase  interval  is  similar  in  all  cases  where  Pe>  - 
22°  and  commonly  ranges  from  (j)  -  -60°  to  -30°  in  which  a  rise  in  the  noise  peak  is 
associated  with  the  interaction-noise  generation.  This  implies  that,  for  >  -22°  and  -15° 
<  Oj  <  -10°,  the  interaction  period  is  irrelevant  to  changes  in  operating  conditions,  or, 
precisely,  to  changes  of  parameters  in  the  freestream  speed,  the  rotor 's  RPM  and  the 
stator's  Ao A. 

Secondly,  from  the  observation  of  the  noise-peaks’  alignment  for  cases  where  fe  > 
-22°,  the  effective  rotor’s  AoA  fe  is  irrelevant  to  the  timing  of  the  noise  peaks. 
Nevertheless,  the  peak-noise  magnitude  does  increase  as  fie  increases. 

Next,  by  inspecting  the  peak  values  during  the  interaction  period  from  cases  where 
Pe  >  -22°,  it  turns  out  that  there  exists  a  direct  relationship  between  the  peak  value  of  the 
interaction  noise  and  the  local  resultant  velocity  Ut  which  is  defined,  in  Fig.2-3,  by  the 
freestream  speed  Uo  and  the  rotor’s  speed  Ur  at  55%  R  of  the  rotor.  It  is  shown  in  Table 
3-2  and  Fig.  3-6  that  the  peak  value  is  proportional  to  the  square  value  of  U,,  i.e.,  U^. 
This  indicates  that  the  peak  value  shown  in  each  noise  signature  is  directly  proportional 
to  the  local  dynamic  pressure,  q, ,  which  is  equal  to  \  poU}  with  a  constant  po  in  this 
study.  However,  Fujita  and  Kovasznay  (1974)  pointed  out  that  the  peak  value  due  to  the 
rotor-stator  interaction  is  proportional  to  the  freestream  dynamic  pressure  qo  defined  as 
\poUo,  as  shown  in  Eq.(1.3).  This  is  because,  in  Fujita  and  Kovasznay’s  study,  the  only 
freestream  speed  Uo  adopted  was  38  m/s  which  is  close  to  the  value  of  41  m/s  for  the  local 
resultant  speed,  U,  with  a  fixed  rotor’s  speed  at  15.7  m/s  (a  circular-cylinder  rod  rotates 
at  300  rpm).  Since  the  values  of  both  Uo  and  Ut  are  too  close  to  be  distinguished,  it  was 
misleading  to  use  qo  in  their  study.  Therefore,  a  correction  for  both  the  calculation  of  the 
noise  prediction  (Fujita  and  Kovasznay  1974)  and  non-dimensionalization  of  the  measured 
noise  (Ho  and  Kovasznay  1976)  should  be  taken  by  the  local  dynamic  pressure  qt  instead 
of  the  freestream  dynamic  pressure  qo- 
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CHAPTER  4 

STATOR  RESPONSE  TO  FLOW  DISTURBANCES 


As  mentioned  in  Chapter  1,  Fujita  and  Kovasznay  (1974)  studied  the  rotor-stator 
interaction  noise  in  a  free  space.  They  simplified  Curie’s  supplementary  formula  to  relate 
a  transient  stator  response  (i.e.,  the  lift  fluctuations)  to  the  corresponding  free-space, 
near-field  dipole  sound  (Eq.  (1.3)).  In  this  study,  the  sound  field  were  measured  in  a 
confined  region  by  the  pipe.  Accordingly,  the  Fujita-Kovasznay  formula  (Eq.(1.3))  may 
not  be  directly  applied  in  the  current  circumstance.  Nevertheless,  the  local  flow 
interaction  results  in  a  net-force  variation  acting  on  the  stator,  which  induces  the  noise 
generation.  In  other  words,  the  noise  is  produced  by  the  change  of  both  and  dC^/  dt, 
but  the  proportional  constant  in  the  equation  will  not  hold. 

4.1.  Stator  Response 

Stator  response  within  a  complete  period  based  on  BPF  is  the  time  trace  of  the 
resultant  aerodynamic  force  on  the  stator  surface.  Since  the  flow  incidence  relative  to  the 
stator  is  varied  with  time,  hence  it  is  hardly  to  determine  the  instantaneous  lift-  and  drag- 
forces’  directions  case  by  case.  Therefore,  the  stator  response  is  resolved  into  the  normal 
and  tangential  components  instead  of  the  lift  and  drag  for  investigation.  In  this  chapter, 
we  are  about  to  review  36  cases  for  the  relationship  between  the  stator  response  and 
interaction  noise  in  a  range  of  -22°  and  -15“  <  a,  <  10°. 

4.1.1.  Operating  Conditions 

In  Section  3.1,  we  have  one  hundred  and  twenty  operating  conditions  to 
investigate  the  interaction-noise  signatures.  In  order  to  examine  the  dependency  between 
both  the  acoustic  fields  and  stator  responses,  thirty-six  conditions  for  the  surface  pressure 
measurements  are  chosen  from  them,  and  are  described  as  follows. 

There  are  six  different  rotor’s  effective  AoA  /?  which  are  -10.30°,  1.32°,  5.82°, 
9.61°,  -21,78°  and  2.20°  with  their  corresponding  incident  velocities  U,  at  55%R  of  9.49, 
13.46,  16,29,  19.92,  11.38  and  20.89  m/s,  respectively,  as  listed  in  Table.  3-1.  Also  from 
Table  3-1,  the  first  four  values  of  are  set  up  by  four  different  rotor’s  RPMs  at  1500, 
2400,  3000  and  3750,  with  a  fixed  freestream  velocity  Ug  at  5.62  m/s,  which  present  four 
corresponding  rotor’s  velocities  at  55%R  of  7.65,  12.23,  15.29  and  19.11  m/s, 
respectively.  However,  the  latter  two  values  of  /?,,  -21.78°  and  2.20°,  correspond  to  the 
two  rotor  RPMs  at  1500  and  3750,  in  turn,  with  a  fixed  freestream  velocity  Ug  at  8.43 
m/s.  Based  on  these  values,  we  will  survey  the  relationship  between  the  stator  response 
and  interaction  noise  for  the  cases  where  -22°  <P^<  10°. 

With  six  stator’s  AoA,  a,,  at  -15°,  ±10°,  ±5°  and  0°,  we  hence  have  36  operating 
conditions  for  a  comprehensive  investigation. 
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4.1.2.  Force  Coefficients  :  Q  &  Q 


To  evaluate  the  components  of  the  stator  response,  the  stator’s  pressure 
distributions  on  both  sides  must  be  integrated.  It  is  straightforward  to  choose  the  positive 
direction  for  the  tangential  force  component  as  pointing  toward  the  trailing  edge  in 
parallel  with  the  stator’s  chord  wise  direction.  Since  the  stronger  impact  happens  on  the 
pressure  side  during  the  flow  interaction  for  the  cases  that  -22°  <  /?  <  0°,  the  normal  force 
component  is  calculated  by  subtracting  the  resultant  force  on  the  suction  side  from  the 
one  on  the  pressure  side;  hence  its  positive  direction  shall  be  normal  to  the  chordwise 
direction  and  pointing  outward  from  the  pressure  side.  A  schematic  description  for  the 
senses  of  both  force  components  is  illustrated  in  Fig.  4-1,  where  N  and  T  represent  the 
normal  and  tangential  forces,  respectively. 

From  Chapter  3,  we  know  that  the  peak-noise  value  is  proportional  to  the  local 
dynamic  pressure  q,  for  the  entire  surveying  range.  We  may  define  the  normal  force 
coefficient,  Q,  and  the  tangential  force  coefficient,  C^,  as  follows: 


C 


N 


N 

q,A 


N 

-c-l 


(4.1a) 


and 


T 

hpU^-c.\ 


(4.1b) 


where  N  and  T  are  magnitudes  of  the  normal  and  tangential  forces,  respectively,  q,  is 
defined  as  with  the  local  resultant  velocity  U,,  and  A  is  the  planform  area.  Since, 

at  present  study,  the  stator  response  is  calculated  on  the  blade  section  at  55%R  of  the 
rotor,  the  planform  area.  A,  for  a  certain  section  of  the  stator  blade  is  given  by  a  product 
of  its  chord,  c,  and  a  unit  span  in  Eq.(4.1). 

Accordingly,  the  normal  and  tangential  forces,  and  acting  on  the  stator  are 
calculated  and  plotted  for  thirty  six  cases  in  Fig.  4-2  through  4-7.  We  shall  discuss  these 
results  in  two  categories:  -22°  <  <  0°  and  0°  <  10°. 


Fluctuations  of  and  Cj.  curves  in  -22"  <  p^<(f 

There  exists  a  resemblance  among  the  curves  in  the  category  of  -22°  <  /?  <  0°, 
as  shown  in  Fig.  4-2  and  4-3.  In  general,  a  rising  slope  in  each  Q  curve  is  consistently 
observed  during  the  interaction  period  from  <f)  =  -60°  to  -30°,  and  ends  up  with  reaching  a 
maximum.  However,  it  is  hardly  to  find,  from  these  figures,  a  unique  trend  from  the 
variations  of  the  Q  curves  in  this  category. 

Meanwhile,  from  Fig.  4-2  &  4-3,  we  observe  that  the  fluctuation  magnitude  of 
each  Q  curve  is  greater  than  that  of  the  corresponding  Q  curve.  This  shows  that  the 
majority  of  surface  pressure  fluctuation  contributes  to  the  normal  force  rather  than  to  the 
tangential  one.  In  other  words,  the  sound  radiation  from  the  stator  surface  is  mainly 
consequent  on  the  normal  force  fluctuation.  The  discussion  will  be  subsequently  focused 
on  the  normal  force  coefficient  Q  for  this  category. 


Fluctuations  of  and  C^.  curves  in  ff 

Both  the  Q  and  Cj.  curves  in  the  category  of  tf  <p^<  10°  are  plotted  from  Fig,  4- 
4  to  4-7.  For  the  smaller  cases,  i.e.,  P^  =  1.32°  and  2.20°,  a  similar  tendency  of  rising 
slopes  in  the  Q  curves  occurs  during  the  interaction  period  in  Fig.  4-4  &  4-5.  It  is 
noticed  that  the  magnitudes  of  these  rising  slopes  are  apparently  smaller  than  those  for  - 
22°  <  /?  <  0°  in  Fig.  4-2  &  4-3.  However,  this  tendency  in  the  curves  is  invalid  for  the 
higher  p^  cases.  For  p^  =  5.82°,  the  rising  slopes  in  the  Q  curves  become  vague  during 
the  interaction  period  in  Fig.  4-6,  and  turn  themselves  into  declining  slopes  for  p^  =  9.61° 
in  Fig.  4-5.  Besides,  the  magnitudes  of  the  curves  for  these  two  cases  are  smaller  than 
those  in  the  smaller  p^  cases  in  Fig.  4-4  &  4-5. 

For  the  curves  in  this  category,  there  does  not  exist  a  consistent  trend  from 
their  variations  on  Fig.  4-4  through  4-7.  Meanwhile,  from  these  figures,  we  observe 
again  that  the  fluctuation  magnitude  of  each  Q  curve  is  greater  than  that  of  the 
corresponding  Cj.  curve,  which  means  that  the  normal  force  indicates  the  main  part  of  the 
stator  response. 

The  above  observations  imply  that  the  stator  response  in  the  category  of  0°  <  /?,  < 
10°  is  weak  and  not  responsible  for  the  majority  of  the  interaction-noise  generation, 
especially  in  the  higher  p^  cases. 


4.1.3.  Fluctuations  :  Cf.  & - 

a^  dt 

Since  sound  radiation  results  from  the  pressure  fluctuation,  a  new  coefficient, 
Cfj ,  will  be  introduced  in  the  following  sections  to  represent  the  fluctuation  part  of  the 
normal  force  coefficient  Q .  Accordingly,  the  time  traces  of  the  new  coefficient,  ,  are 
plotted  in  part  (a)  of  Fig.  4-8  to  4-13.  Each  Q  curve  is  simply  the  fluctuation  portion  of 

its  corresponding  coefficient  Q,  except  that  the  Q  fluctuations  for  p^  =  9.61°  in  Fig.  4-7 
are  inverted  and  plotted  in  Fig.  4-13.  This  exception  is  made  since  we  have  no 
knowledge  about  the  influence  of  the  force  orientation  to  the  directivity  of  the  acoustic 
field  (note  that  we  only  know  that  the  acoustic  wave  arriving  at  0.7  m  downstream 
becomes  a  plane  wave.)  Hence,  in  Fig.  4-13,  the  inverted  Cf^  curves  show  a  consistent 
rising  slope  during  the  interaction  period  from  ^  =  -60°  to  -12°. 

However,  the  peaks  of  Cf^  in  part  (a)  of  Fig.  4-8  to  4-13  do  not  coincide  with  the 
corresponding  peak  noises  observed  in  Fig.  3-2  (b),  (c),  (d)  8c  (e)  and  Fig.  3-3  (b)  &  (e). 
We  will  discuss  it  later  in  this  section.  _ 

Another  fluctuation  item  in  which  we  are  interested  is  the  time  derivative  of 
itself.  From  Fujita-Kovasznay  simplified  formula  in  Chapter  1,  the  lift- fluctuation  term 
contributes  mainly  to  the  near-field  noise,  while  its  time-derivative  term  dominates  the 
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far- field  noise  (Lighthill,  1962).  Meanwhile,  since  we  have  the  normal  force  instead  of 
the  lift  force,  the  Fujita-Kovasznay  formula  is  rewritten  in  terms  of  Q  as  following: 


S.P.  =  k,{  c„  +  }  , 

Oq  at 
near  field  far  field 


(4.2) 


where  .  cos6 ,  as  defined  in  Chapter  1.  Note  that,  given  /?  =  0.7  m  and  =  340 

AtiR  ° 


m/s,  then  R  !  h  0.00206  sec.  Time  traces  of 
plotted  in  part  (b)  of  Fig.  4-8  through  4-13. 


R  dC^ 

«o  dt 


are  subsequently  obtained  and 


First,  from  these  figures,  both  fluctuation  magnitudes  of  the  Q  and 


RdC^ 

Uq  dt 


curves  decrease  as  decreases.  Secondly,  for  a  fixed  the  magnitude  of  Q  is  of  the 


same  order  as  that  of 


dt 


■.  Therefore,  the  term  in  Eq.(4.2)  contributes  as  much 


to  the  sound  radiation  from  the  stator  surface  as  the 


R  dC^ 

flo  dt 


Thirdly,  we  also  observe  that  the  rising  slope  occurring  in  each 


term  calculated  at  0.7  m. 
R  dC. 


«o  dt 


curve  during 


the  interaction  period  is  more  correlated  with  the  peak  noise  from  the  measurements  than 
that  in  the  corresponding  curve.  This  observation  is  in  agreement  with  the  fact  that 


R  dC^ 

«o  dt 


governs  the  far-field  noise  for  the  free  space  as  R  approaches  infinity,  as 


described  in  Chapter  1 . 

In  the  light  of  Eq.(4.2),  we  can  calculate  the  resultant  contribution  from  both 
fluctuation  terms  which  will  be  discussed  in  the  next  section. 


4.2.  Source  Term  vs.  Measured  Noise 

We  define  the  source  term  by  the  summation  of  both  terms  in  curly  brackets  on 

R  dC^ 

the  nght-hand  side  of  Eq.(4.2),  i.e.,  { -i- - - — }.  Then,  the  source-term  curves  are 

Oq  at 

plotted  with  the  measured  noise  signatures  from  Fig.  4-14  to  4-19. 
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For  -22°  <  P^<Q°  in  Fig.  4-14  &  4-15,  there  is  a  good  agreement  in  peak  value 
between  both  the  source-term  and  measured  results  during  the  interaction  period.  It  is 
surprising  because  the  sound  propagation  in  the  pipe  is  not  influenced  so  much  due  to  the 
confined  geometry  as  we  expected.  For  (f  <p^<  10°  in  Fig.  4-16  to  4-19,  the  peak  values 
in  the  source-term  results  become  smaller  than  those  in  the  measured-noise  results 
specifically  as  increases. 

Note  that  the  source  term  is  calculated  based  on  the  stator  response.  Therefore, 
the  stator  response  becomes  the  dominant  resource  for  the  interaction-noise  generation  if 
the  rotor’s  effective  AoA,  is  negative.  If  p^  is  positive,  it  becomes  less  responsible  for 
the  majority  of  the  interaction-noise  generation  with  increasing  p^. 


4.2.1.  Amplitude  factor  k  v^.  interaction  angle  y 

An  amplitude  factor,  k,  is  defined  as  a  ratio  of  the  source-term  peak  value  to  the 
peak  noise,  and  can  be  easily  obtained  by  measuring  and  comparing  the  peak  values  from 
Fig.  4-14  to  4-19,  during  the  interaction-noise  generation.  For  each  k,  there  exists  both 
the  corresponding  values  of  and  /?  which  enable  us  to  define  an  interaction  angle,  y, 
via  a  simple  triangular  relationship  illustrated  in  Fig.  4-20. 

In  Chapter  3,  we  found  that,  from  Eq.(3.4),  the  noise-peak  value  is  also 
proportional  to  the  freestream  velocity,  U^,  for  the  cases  that  >  -22°.  In  addition,  as 
mentioned  in  Section  6.1,  the  thirty-six  operating  conditions  for  surface  pressure 
measurements  are  chosen  based  on  two  freestream  velocities  at  =  5.62  and  8.43  m/s. 
The  diagrams  are,  hence,  plotted  for  both  of  them  in  Fig.4-21(a)  &  (b),  respectively. 

Generally  in  Fig.4-21(a)  &  (b),  for  a  fixed  stator’s  AoA  a^,  the  amplitude  factor  k 
increases  as  the  interaction  angle  y  increases.  According  to  the  definition  of  k,  the  higher 
the  K,  the  stronger  the  noise  generation  from  the  stator  response.  From  Fig.4-20,  for  a 
fixed  a,,  the  rotor’s  effective  AoA  p^  decreases  while  the  interaction  angle  increases. 

Note  that,  as  x:  approaches  unity  for  -22“  <  p^<0°  in  Fig.  4-21  (a)  &  (b),  there 
exists  a  high  dependency  between  both  the  peak  values  of  the  source  terms  and  measured 
noises,  which  suggests  that  the  stator  response  is  the  major  resource  to  the  noise 
generation.  Since  the  source  term  is  calculated  from  the  pressure  distributions  on  the 
stator’s  blade  section  at  55%R,  the  good  agreement  between  both  of  the  peak  values 
implies  that  the  stator  response  at  55%R  is  typical  of  the  entire  stator  response  along  the 
span  with  a  factor  k,  essentially. 

For  0“  <  /?,  <  10°,  the  value  of  x- decreases  as  p^  increases  with  a  fixed  a,  in  Fig.  4- 
21  (a)  &  (b).  In  other  words,  the  stator  response  in  this  category  does  not  contribute  so 
much  to  the  noise  generation  as  that  in  -22“  <  P,<  0“.  Hence,  it  implies  that  a  transient 
rotor  response  becomes  responsible  for  the  majority  of  interaction-noise  generation  for  0“ 
<P<\Q°. 


28 


4.3.  Summary 

In  this  chapter,  we  may  summarize  as  follows. 

•  In  the  light  of  the  Fujita-Kovasznay  formula  [Eq.(4.2)],  we  find  that  the  source  terms 

calculated  from  the  stator’s  surface  pressure  distributions  agree  with  the  measured 
noise  signatures  for  the  category  of  -22”  <  /?  <  0”,  while,  for  0°  <  <  10°,  this 

dependency  becomes  less  as  /?,  increases. 

•  Meanwhile,  by  defining  two  parameters,  i.e.,  the  amplitude  factor  ic  and  the 
interaction  angle  y,  we  have  reached  the  following  conclusion;  For  a  fixed  stator’s 
AoA,  the  interaction-noise  generation  highly  depends  upon  the  stator  response  as 
the  interaction  angle  y  increases,  while,  if  y  decreases,  the  rotor  response  shall  be 
mainly  responsible  for  the  interaction-noise  generation. 


CHAPTER  5 

FLOW  FIELD  AROUND  THE  STATOR 


One  of  the  main  tasks  in  the  present  study  is  to  understand  the  physical 
mechanism  of  the  rotor-stator  interaction-noise  generation  by  correlating  both  the 
acoustic  field  (ejfect)  and  the  corresponding  aerodynamic  field  (cause).  In  this  chapter, 
we  will  examine  the  local  flow  fields  in  detail  and  qualitatively  correlate  them  to  the 
interaction  noise.  Note  that  the  active  interaction  phase  interval  for  the  cases  where  >  - 
22°  is  previously  identified  as  the  phase  angle  ^  from  -60°  to  -30°.  Nevertheless,  we  may 
as  well  extend  the  investigation  interval  from  <j)  =  -60°  up  to  -12°  for  study  in  order  to 
understand  more  details  about  the  unsteadiness  of  the  interaction  flow  field. 

Meanwhile,  since  the  peak  value  of  the  interaction  noise  increases  with  the  rotor’s 
effective  AoA  for  cases  where  >  -22°  (Chapter  3),  it  is  convenient  to  investigate  the 
flow  field  categorized  by  p^  in  this  chapter. 

5.1.  Operating  Conditions  for  Surface  Pressure  and  Velocity  Measurements 

There  are  thirty-six  operating  conditions  for  the  surface  pressure  measurements 
which  are  same  as  described  in  sub-section  4.1.1. 

For  the  velocity  measurements,  sixteen  conditions  are  selected  from  those 
conditions  in  surface  pressure  measurements.  Except  only  choosing  =  0°  for  case  of  P^ 
=  -21.78°,  three  a/s  are  set  up  at  0°  and  ±10°  for  the  remaining  five  P/s,  i.e.  -10.30°, 
1.32°,  5.82°,  9.61°  and  2.20°.  Fig.  5-1  shows  three  meshes  of  hot-wire  stations  at  a,  = 
10°,  0°  and  -10°  used  for  velocity  measurements  around  the  stator’s  FIRST  30%C.  Each 
mesh  has  approximately  4,500  hot-wire  stations. 

In  the  following  sections,  the  results  of  these  two  measurements  will  be  discussed 
and  correlated  to  the  noise  results  in  order  to  survey  the  noise-generation  mechanism  due 
to  rotor-stator  interaction  for  cases  where  -22°  <P^<  10°. 


5.2.  Categorization  of  Interaction  Flow  Fields  by 

Recall  that,  in  sub-section  3.2.3,  two  typical  noise  patterns  have  been  observed 
within  -22°  <  P/<  10°  and  distinguished  by  cases:  p^  =  ±10°.  Moreover,  from  sub-section 
3.2.1,  it  has  been  pointed  out  that,  for  cases  where  0°  <  p^<  10°,  the  noise-peak  value 
increases  as  the  stator’s  AoA  a,  decreases.  Therefore,  the  zero  degrees  of  p^  which 
represent  the  rotor’s  zero-lift  AoA  will  be  a  somewhat  rough  but  good  dividing  line  for 
categorizing  the  cases.  Accordingly,  we  may  easily  survey  the  interaction  flow  fields 
with  two  categories:  -22°  <  /?^  <  0°  and  0°  <  <  10°. 
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In  the  following,  two  cases,  i.e.  >0^  =  -10.30°  and  9.61°,  will  be  utilized  to 
represent  both  categories  for  investigation.  The  experimental  results  regarding  the  rest 
cases  can  be  found  in  Hu  and  Ho’s  report  (1997). 

5.3.  The  Flow  Fields  for  Cases:  -22°  <  0°. 

The  cases  for  investigation  in  -22°  <  <  0°  are  =  -10.30°  with  -15°  <  a  <  10°. 

5.3.1.  Characteristics  of  surface  pressure  results  for  =  -10.30° 

Unsteady  surface  pressure  distributions  for  =  -10.30°  are  shown  in  Fig.  5-2, 
with  six  stator’s  AoA  during  -60°  <  <j)<  -12°.  The  pressure  distributions  in  Fig.  5-2  are 
non-dimensionalized  by  the  local  dynamic  pressure  q,,  i.e.  and  plotted  against 

the  normalized  chordwise  coordinate,  x/c.  As  mentioned  previously,  the  pressure 
distributions  on  both  sides  are  shown  up  to  85%C  due  to  the  limitation  of  measurements 
(Chapter2). 

The  physical  significance  of  Fig.5-2  is  described  as  follows.  First,  for  all  six 
stator’s  AoA  a/s,  the  pressure  distributions  on  the  stator’s  pressure  side  fluctuate  more 
drastically  than  those  on  its  suction  side  during  the  interaction  period.  This  seems 
consequent  to  the  fact  that  the  pressure  side  is  the  side  facing  to  the  rotor’s  approach 
(Fig.  2-3)  and  is  supposed  to  experience  the  most  impact  during  the  interaction  period. 

Secondly,  we  will  look  into  how  the  surface  pressure  varies  during  the  interaction 
period  on  the  stator’s  pressure  side.  For  a,  >  0°,  which  means  the  stator  is  tilted  to  its 
pressure  side  for  a  positive  angle  of  attack,  most  of  the  pressure  side  experiences  negative 
pressure,  shown  in  Fig.  5-2(a)  &  (b),  and  the  dramatic  change  observed  is  a  diminishing 
region  of  negative  pressure  gradient  from  20%C  down  to  2%C  due  to  a  decrease  of  the 
surface  pressure.  For  a,  =  0°,  a  similar  progress  is  observed  in  Fig.  5-2(c).  The  difference 
is  that  the  negative-pressure-gradient  region  shrinks  from  45%C  at  ^  =  -60°  down  to 
5%C  at  ^  =  -12°,  which  causes  a  strong  pressure  drop  on  the  pressure  side.  For  <  0°, 
the  pressure  distributions  on  the  pressure  side  are  mostly  positive,  as  shown  in  Fig.  5- 
2(d),  (e)  &  (f),  and  the  highest  pressure  indicates  the  stator’s  front  stagnation  point  is 
observed  and  located  within  the  first  5%C  of  the  pressure  side.  In  a  like  manner,  a 
diminishing  and  decreasing  negative-pressure-gradient  region  occurs  from  50%C  down 
to  10%C  on  this  side  within  a  phase  interval:  (f>  =  -60°  to  -12°. 

These  results  show  that,  for  all  six  stator’s  AoA  a/s  with  =  -10.30°,  a 
negative-pressure-gradient  region  on  the  stator’s  pressure  side  diminishes  during  the 
rotor-stator  interaction  and  is  identified  as  the  active  interaction  region  on  the  stator. 

5.3.2  The  flow  fields  for  p^  =  -10.30°  at  six  a/s 

It  is  noticed  that  this  pressure  decrease  is  simultaneously  accompanied  by  a 
sudden  drop  of  the  negative  pressure  gradient  in  the  same  active  region  for  the  case:  P^  = 
-10.30°  with  all  six  a/s.  This  active  region  generally  ranges  from  the  stator’s  leading 
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edge  up  to  30%C  on  its  pressure  side.  Hence,  the  following  velocity  results  will  be 
confined  in  this  active  region  (Fig.  5-1).  We  should  inspect  the  cause  of  the  drastic  drop 
in  both  the  surface  pressure  and  its  corresponding  negative  gradient  during  the  interaction 
period.  This  discussion  will  be  initiated  with  the  case;  =  10°  for  =  -10.30°,  and 
extended  to  the  rest  a’s  at  5“,  0°,  -5“,  -10“  and  -15“. 

5.3.2.1  Case  :  =  10°. 

Fig.  5-3  shows  a  history  of  velocity  contours  for  a  =  10°  and  p^  =  -10.30°  from 
phase  angle  ^  =  -60°  to  -12°.  Each  velocity  contour  at  a  certain  phase  angle  is  cross- 
plotted  from  phase-averaged  time-trace  data  measured  at  different  hot-wire  stations  (Fig. 
5- 1(a)).  It  should  be  borne  in  mind  that,  since  the  velocity  measurements  were  done  by  a 
single  hot-wire  probe  which  measures  the  heat-transfer  effect  in  the  normal  direction  to 
the  hot  wire,  hence  the  magnitude  U(x,y,^)  of  velocities  measured,  as  shown  in  Fig.  5-3, 
is  the  resultant,  (u  +  v)"^,  of  two  velocity  components,  u  and  v,  in  a  Cartesian  or  x-y 
coordinate  system  as  plotted  in  Fig.  5-3.  This  indicates  that  the  contour  lines  shown  in 
Fig.  5-3  are  not  the  instantaneous  streamlines  so  that  we  can  not  apply  them  to  identify 
the  local  flow  directions.  Despite  this  drawback,  these  velocity-contour  plots  still  provide 
us  profound  information  about  the  interaction  flow  field  as  pictures  obtained  by  flow- 
visualization  methods. 

We  shall  investigate  the  velocity-contour  plots  in  Fig.  5-3  in  parallel  with  the 
corresponding  surface  pressure  distributions  from  Fig.  5-2(a).  At  ^  =  -60°  in  Fig.  5-3,  the 
rotor’s  wake  which  is  shaded  by  a  low-speed  region  is  clearly  entering  in  the  left-bottom 
comer  of  the  velocity-contour  plot  in  Fig.  5-3.  We  also  observe  that,  at  this  moment, 
another  low-speed  region  occurs  around  the  leading  edge  on  the  stator’s  suction  side. 
This  low-speed  region  is  recognized  as  the  instantaneous  front-stagnation- flow  regime  on 
the  stator,  and  is  also  validated  by  the  existence  of  a  maximum  and  positive  pressure  on 
the  suction  side  at  ^  =  -60°  in  Fig.  5-2(a).  This  front  stagnation  point  located  on  the 
suction  side  is  subsequently  on  the  positive  stator’s  AoA  a,  which  is  rotated  10°  to  the 
pressure  side.  On  the  pressure  side,  there  exists  a  high-speed  spot  marked  by  a  10  m/s 
contour  line  from  the  leading  edge  to  approximate  10%C.  As  we  will  see  later,  this  high¬ 
speed  spot  is  a  strongly  accelerated  flow  region  which  is  an  inviscid  effect  caused  by  the 
separated  rotor’s  wake  approaching  to  the  stator  in  the  present  case.  In  this  flow 
acceleration  region,  a  pressure  drop,  i.e.  a  negative  pressure  gradient  region,  is 
consistently  observed  at  ^  =  -60°  on  the  pressure  side  in  Fig.  5-2(a). 

Next,  at  ^  =  -54°,  the  flow  field  does  not  change  too  much  by  observing  both  Fig. 
5-2(a)  and  5-3,  except  that  the  rotor’s  wake  moves  upwards  and  appears  clearer.  The 
wake  behind  the  rotor  is  separated  due  to  a  negative  rotor’s  AoA  p^  of  -10.30°.  At  this 
circumstance,  the  rotor  itself  serves  as  a  bluff  body,  and  the  boundary  layer  on  its 
pressure  side  has  already  separated  to  form  a  thick  wake  behind. 

Between  ^  =  -48°  and  -12°,  it  is  observed  that  a  strong  interaction  happens  in  the 
flow  field  by  inspecting  the  consecutive  velocity  contours  in  Fig.  5-3.  At  ^  =  -48°,  the 
inviscid  wake  behind  the  rotor  apparently  approaches  to  the  stator’s  pressure  side.  The 
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contour  lines  become  squeezed  within  a  narrowing  region  formed  by  the  margin  between 
the  rotor’s  wake  and  the  stator’s  pressure-side  surface,  especially  near  the  stator’s  leading 
edge.  This  indicates  that  around  the  highly-curved  leading  edge  portion  there  occurs  a 
stronger  velocity  gradient  at  the  present  phase  than  that  at  the  previous  phase. 
Meanwhile,  it  is  also  noticed  that  the  10  m/s-high-speed  spot  on  the  pressure  side  shrinks 
at  ^  =  -48°  due  to  the  rotor’s  wake  approaching. 

Three  sequential  velocity-contour  plots  af  (j>  =  -42°,  -36°  and  -30°  clearly  show 
that  an  interaction  process  takes  place  as  the  rotor’s  wake  (the  separated,  inviscid  flow) 
breaks  in  the  stator’s  boundary  layer  (an  attached,  viscous  flow).  This  process  is  the 
impinging  wake  phenomenon  identified  by  Fujita  and  Kovasznay  (1974).  During  the 
interaction  process,  we  observe  that  the  rotor’s  wake  starts  to  disappear  and  to  connect 
itself  to  the  stator’s  leading-edge  flow  regime.  This  suggests  that  an  active  region  on  the 
stator’s  surface  will  erupt  from  its  originally  attached  boundary  layer  due  to  the 
impingement  of  the  rotor’s  wake. 

Boundary  layer  at  stator’s  leading  edge 

This  eruption  of  the  laminar  boundary  layer  is  validated  by  an  inspection  of  time 
variations  of  velocity  distributions  along  two  lines:  one  line  is  normal  to  the  leading  edge 
and  the  other  normal  to  the  point  at  2%C  on  the  pressure  side,  as  illustrated  in  Fig.  5- 
4(a).  The  temporal  variations  of  velocity  contours  along  the  specific  line  normal  to  the 
leading  edge  is  plotted  in  Fig.  5-5(a),  while  the  normal  coordinate,  n,  has  been  non- 
dimensionalized  by  the  gap  distance  of  0.71cm  at  55%R.  Meanwhile,  Fig.  5-6(b) 
shows  a  portion  of  Fig.  5-6(a)  which  is  10%  of  d^^  from  the  leading  edge.  The  velocity- 
contour  values  are  also  non-dimensionalized  by  the  freestream  velocity  ,  5.62  m/s  for 
the  present  case,  in  Fig.  5-6(a)  and  (b),  where  the  solid  contour  lines  represent  the  values 
greater  than  or  equal  to  one  and  the  dashed  lines  values  less  than  one. 

Since  both  the  rotor’s  and  stator’s  positions  are  perpendicular  to  each  other  at 
phase  angle  ^  =  +180°,  the  flow  field  around  the  stator’s  leading  edge  portion  in  Fig.  5-4 
(b)  looks  like  the  stator  sitting  in  a  uniform  stream  with  less  disturbance  from  the  most 
faraway  rotor’s  position  within  a  revolution  cycle.  Hence,  before  the  strong  interaction 
happens  at  about  (j>  =  -60°,  the  contour  line  with  value  one  in  Fig.  5-5(b)  indicates  the 
boundary  layer  thickness  at  the  stator’s  leading  edge  in  a  certain  sense,  and  has  an 
estimated  thickness  with  an  edge  velocity  U„?it  (p  =  -180°  is  about  0.015  d^^,  as  observed 
in  Fig.  5-5(b). 

From  ^  =  -180°  to  -60°  in  Fig.  5-5(b),  the  boundary  layer  at  the  stator’s  leading 
edge  is  attached  and  laminar.  It  is  noticed  that  this  boundary  layer  is  suppressed  by  the 
presence  of  a  flow  acceleration  process  on  the  leading  edge  marked  by  a  set  of  centered 
solid  lines  from  ^  =  -105°  to  -60°,  in  Fig.  5-5  (a)  &  (b).  The  flow-acceleration  process 
ends  with  the  boundary-layer  thickness  reaching  a  minimum  value  of  0.01  d^^  aX  (/>  =  -60° 
in  Fig.  5-5(b).  This  is  also  supported  by  an  observation  of  a  10  m/s-high-speed  spot  from 
the  leading  edge  to  10%C  on  the  stator’s  pressure  side  at  ^  =  -60°  in  Fig.  5-3.  This  flow- 
acceleration  process  is  due  to  a  contraction  formed  instantaneously  by  the  rotor’s 
separated  wake  margin  and  the  curved  leading  edge  geometry  on  the  stator’s  pressure 
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side  (as  illustrated  in  Fig.  5-7(a)).  A  flow  deceleration  on  the  leading  edge  is  observed 
from  (j)  =  -60°  to  -42°  in  Fig.  5-5(b),  while  the  boundary  layer  grows  a  little  during  this 
process.  This  deceleration  process  is  also  validated  with  the  diminishing  of  the  10  m/s- 
high-speed  spot  from  (j)  =  -60°  to  -42°  in  Fig.  5-3.  After  ^  =  -42°,  the  boundary  layer  at 
the  stator’s  leading  edge  suddenly  erupts  and  extends  along  the  normal  direction  to  the 
rotor’s  surface,  as  shown  in  Fig.  5-5(a)  &  (b).  This  boundary  layer  eruption  is  due  to  an 
impingement  of  the  rotor’s  wake  onto  the  stator’s  surface,  in  Fig.  5-3  and  Fig.  5-5(a)  & 
(b).  This  process  is  a  viscous-inviscid  interaction  between  a  viscous,  attached  flow  (the 
stator’s  boundary  layer)  and  an  inviscid,  separated  flow  (the  rotor’s  separated  wake) 
(Didden  and  Ho  1985). 

In  Fig.  5-2(a),  we  observe  that  a  pressure  drop  is  accompanied  with  a  decrease  of 
negative  pressure  gradient  (i.e.,  dpjdsKQ)  at  first  5%C  of  the  stator’s  pressure  side. 
Since  a  pressure  gradient  along  the  surface  is  necessary  to  sustain  a  flux  of  vorticity  into 
the  fluid  (Panton  1984),  the  Navier-Stokes  equation  along  the  streamwise  direction  on  a 
solid  surface  for  a  2-D  flow  becomes 


5©^ 

dn 


|i.  ds 


on  the  surface 
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where  is  the  vorticity  gradient.  Hence,  both  time  traces  of  the  normalized  surface 


pressure  and  the  normalized  surface  vorticity  gradient,  dw*/dn*  (calculated  by 
-\l\i{dplds)-{dg^py jU,),  at  the  leading  edge  are  plotted  in  Fig.  5-5  (c)&(d), 
respectively. 

Fig.  5-5(c)  and  (d)  show  that,  from  (j)  =  -60°  to  -30°,  a  sudden  pressure  drop  at  the 
leading  edge  (Fig.  5-5(c))  obviously  accompanies  with  a  vorticity  flux  (i.e.,  an  increase 
of  vorticity  gradient)  into  the  fluid  flow  (Fig.  5-5(d)).  Accordingly,  a  strong  flow 
interaction  between  the  rotor  and  stator  happens  at  ^  =  -60°,  from  which  proceeds  a 
sudden  flux  of  vorticity  into  the  fluid  around  the  stator’s  leading  edge  portion. 

Moreover,  the  boundary  layer  eruption  at  the  leading  edge  continues  until  the 
phase  angle  ^  =  0°  after  which  the  boundary  layer  is  gradually  recovered  by  the  rotor’s 
passing,  and  it  reattaches  to  the  stator’s  leading  edge  after  (j)  =  30°  in  Fig.  5-5(a)  &  (b). 


Boundary  layer  at  2%C  on  stator’s  pressure  side 

On  a  similar  basis.  Fig.  5-6(a)  &  (b)  shows  the  time-trace  contour  along  the 
normal  line  originated  from  the  point  at  2%C  on  the  stator’s  pressure  side.  Since  the  flow 
is  accelerated  by  the  favorable  pressure  gradient  to  overcome  the  curved  geometry 
around  the  leading  edge  at  =  -180°  in  Fig.  5-4(b),  an  edge  velocity  of  the  boundary 
layer  at  the  current  point  should  be  larger  than  that  at  the  leading  edge.  It  is  such  a  short 
distance  from  the  leading  edge  to  the  current  point  at  2%C  that  the  laminar  boundary 
layer  thickness  over  this  distance  should  grow  very  little  at  ^  =  -180°.  Hence,  at  ^  =  - 
180°  in  Fig.  5-6(b),  the  contour  line  with  a  constant  value  of  1.5  is  capable  of  indicating 
a  boundary  layer  thickness  of  0.015  d^^^  (note  that  an  approximate  value  of  Ug  at  the  edge 
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of  the  leading-edge  boundary  layer  has  been  adopted  earlier  with  0.015  as  a  boundary 
layer  thickness  at  ^  =  -180°.)  From  Fig.  5-6(b),  we  observe  that  a  thin  boundary  layer  at 
this  point  is  mostly  attached  during  one  cycle  except  a  phase  interval  in  which  an 
eruption  of  the  boundary  layer  happens  due  to  the  flow  interaction.  Thus,  following  the 
contour  line  of  1.5  in  Fig.  5-6(b),  the  boundary-layer  thickness  starts  decreasing  at  ^  =  - 
90°  due  to  the  flow  acceleration  above  and  reaches  a  minimum  thickness  at  ^  =  -66°.  In 
the  meantime,  the  rotor’s  wake  clearly  shows  up  at  ^  =  -90°  on  the  top  of  Fig.  5-6(a),  and 
approaches  to  the  point  at  2%C  on  the  pressure  side.  A  flow-deceleration  process  above 
the  occurs  from  ^  =  -66°  to  -42°  in  Fig.  5-6(b),  which  is  also  observed  in  Fig.  5-3.  The 
boundary  layer  grows  up  during  the  flow-deceleration  process,  and  eventually  erupts, 
after  <}>  =  -42°,  due  to  an  impingement  of  the  rotor’s  wake  at  this  point.  The  eruption  of 
the  boundary  layer  lasts  until  about  ^  =  15°  at  which  flow  reattachment  is  observed  in 
Fig.  5-6(b). 

Fig.  5-6(c)  and  (d)  show  the  time  variations  of  surface  pressure  and  vorticity 
gradient  at  the  current  point.  In  these  figures,  a  dramatic  pressure  drop  and  a  sudden  rise 
of  vorticity  gradient  are  observed  again  from  (j)  =  -60°  to  -30°.  Therefore,  it  follows  that 
the  approaching  wake  is  responsible  for  the  sudden  drop  in  the  surface  pressure  about  the 
stator’s  leading  edge,  which  causes  an  interaction  process  between  a  viscous  flow  (the 
stator’s  boundary  layer)  and  an  in  viscid  flow  (the  rotor’s  wake). 

Since  this  viscous-inviscid-interaction  process  happens  exactly  during  a  sudden 
rise  of  the  interaction  noise  in  Fig.  3-2(b),  it  is  directly  responsible  for  the  interaction- 
noise  generation  in  the  present  case.  The  interaction  process  from  (j>  =  -90°  to  -30°  is 
illustrated  in  Fig.  5-7  (a)&(b). 

Post-interaction-noise  period 

After  the  interaction  process  in  Fig.  5-3,  the  rotor’s  wake  is  completely  cut  off  as 
shown  in  the  velocity-contour  plot  at  ^  =  -24°.  Then,  an  attachment  of  the  rotor’s  wake 
onto  the  stator’s  leading-edge  portion  are  shown  in  contour  plots  at  ^  =  -18°  and  -12°. 

5.3.2.2.  A  brief  summary  for  =  -10.30°  with  six  a/s 

The  velocity-contour  plots  for  p^  =  -10.30°  and  a,’s  at  0°  and  -10°  are  shown  in 
Figs.  5-8  and  5-9,  respectively.  Generally  speaking,  the  interaction-noise  generation  is 
due  to  the  surface-pressure  fluctuations  on  the  stator’s  pressure  side  for  all  six  a/s  with 
P^  =  -10.30°,  and  these  fluctuations  are  induced  by  the  flow  separation  around  the  stator’s 
leading  edge  while  the  rotor’s  wake  approaches  to  the  stator’s  pressure  side.  It  is 
interesting  to  point  out  that  the  surface-pressure  fluctuations  are  different  in  two 
categories  (Fig.  5-2).  For  <  0°,  a  negative-pressure  drop  in  the  active  region  on  the 
stator’s  pressure  side  is  esteemed  as  the  surface-pressure  fluctuation,  while  a  positive- 
pressure  drop,  in  contrast,  is  observed  for  a,  >  0°.  An  accompanying  negative-pressure- 
gradient  decrease  is  also  observed  in  each  a/s  case.  This  is  very  important,  as  we  shall 
see  later,  when  we  try  to  seek  a  unique  mechanism  for  the  interaction-noise  generation 
for  -22”  <p^<  0°. 


35 


5.4.  Vorticity  Flux  vs.  Sound  Radiation  for  -22°  <  /3^<0° 


From  the  above  discussion  of  both  the  surface  pressure  and  velocity  results,  it  is 
obvious  that,  around  the  stator’s  leading  edge,  the  impingement  of  the  rotor’s  wake 
causes  locally,  in  the  main,  a  strong  surface  pressure  fluctuation  which  radiates  sound 
wave.  In  other  words,  the  unsteady  stator  response  induced  by  this  pressure  fluctuation  is 
the  source  of  sound  radiation  for  the  category;  -22°  <  /?^  <  0°  and  -15°  <  a,  <  10°.  We 
also  find  that  whatever  this  surface  pressure  fluctuation  varies  case  by  case,  a  drop  of 
negative  surface  pressure  gradient  is  consistently  observed  in  each  case  during  the 
interaction  period.  This  surface  pressure  gradient  is  directly  associated  to  the  surface 
vorticity  gradient,  as  described  in  the  previous  section.  Therefore,  in  the  following,  a 
unique  description  based  on  the  concept  of  vorticity  flux  is  developed  to  summarize  the 
relationship  of  the  negative  pressure  gradient  to  the  stator  response  and  the  sound 
radiation. 


5.4.1.  Concept  of  vorticity  flux 

The  concept  of  vorticity  flux  is  due  to  Lighthill  (1963).  He  pointed  out  that:  "The 
tangential-vorticity  source  strength  has  a  simple  relation  to  pressure  gradient,  at  least 
for  flow  over  a  stationary  plane  surface.  ....  and  because  at  a  solid  surface  transfer  of 
momentum  by  convection  is  absent,  so  that  transfer  by  diffusion  must  exactly  balance 
transfer  by  pressure  gradient.'  Furthermore,  the  balance  between  both  the  vorticity 
generation  and  transport  is  mainly  the  cause  to  change  the  vorticity  field  in  an  unsteady 
flow  (Reynolds  and  Carr  1985,  Shih  and  Ho  1994).  Thus,  for  -22°  <  <  0°,  the 

separated  flow  behind  the  rotor  induces  the  unsteadiness  of  the  vorticity  field  around  the 
stator’s  leading-edge  surface  to  produce  the  vorticity  into  the  fluid  and  hence  the  sound. 

On  the  basis  of  this  statement,  a  relationship  can  be  obtained  between  vorticity 
flux  (i.e.,  the  tangential-vorticity  source  strength)  and  pressure  gradient  on  a  solid 
surface.  First,  we  derive  this  relationship,  and  then  take  advantage  of  it  to  explore  a 
connection  of  the  surface  vorticity  flux  to  the  noise  generation. 

Vorticity  flux  is  used  to  phrase  the  vorticity  generation  across  a  solid  surface  into 
the  flow.  Consider  a  local  s-n-z  coordinate  system  on  the  stator  surface  as  illustrated  in 
Fig.  5-10,  where  s  and  n  represent  the  streamwise  and  the  normal  coordinates  which  are 
tangential  and  perpendicular  to  the  surface,  respectively,  and  z  the  coordinate  which  is  bi¬ 
normal  to  both  s-  and  n-  coordinates  on  the  surface. 

The  z-component  of  a  positive  vorticity  flux  vector  a  across  a  surface  with  the 
normal  n  is  defined  as  (Pan ton  1984) 
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1  dp 
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(5.1) 


Eq.(5.1)  shows  a  relationship  between  vorticity  flux  and  surface  pressure 
gradient.  Since  the  surface  pressure  measurements  have  been  conducted  to  present  study, 
their  results  enable  us  to  calculate  the  vorticity  flux  from  the  surface  pressure  gradient  via 
Eq.(5.1). 

Note  that  is  equal  to  the  negative  value  of  darjdn.  In  other  words, 
ifdp/ds<0  and  5co^/5«>0,  then  a^<0;  a  negative  vorticity  is  generated  from  the 
surface  to  the  fluid.  Namely,  the  surface  is  a  vorticity  source.  Conversely,  if  dplds>0 
and  d(ojdn<0,  then  a,  >0:  a  negative  vorticity  enters  into  the  surface  from  the  fluid; 
hence  the  surface  a  vorticity  sink. 


5.4.2.  Contour  and  Surface  Plots  of  dto  ^  /dn  :  for  -22°  <  <  0° 

Since  the  most  active  region  during  the  interaction  period  is  typically  located  on 
the  first  25%C  of  the  stator’s  pressure  side  for  -22°  <  0°,  hence  we  may  calculate  the 
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vorticity  flux  for  this  specific  region  from  the  surface-pressure  results  to  see  its  variation 
within  one  BPF  harmonic. 

For  =  -10.30°,  the  surface- vorticity-gradient  contour  and  surface  plots  are 
shown  in  Fig.  5-11  at  six  stator’s  AoA  «/s.  Two  non-dimensionlalized  base  axes  plotted 
in  Fig.  5-11  are  the  temporal  and  spatial  coordinates  in  terms  of  the  phase  angle  (j)  and 
the  normalized  stream  wise  coordinate,  s*,  respectively.  Note  that,  in  this  figure,  the 
normalized  streamwise  coordinate,  s*,  is  confined  from  the  stator’s  leading  edge  up  to 
25  %C  of  the  pressure  side,  which  contains  the  active  region  for  the  flow  interaction.  The 
contour  plot  is  stacked  upon  the  surface  one  for  the  convenience  of  inspecting  the 
variation  of  vorticity  flux  within  a  complete  cycle  in  the  active  region.  The  height  of  both 

contour  and  surface  plots  represents  ddi'Jdn  ,  the  dimensionless  magnitude  of  surface 
vorticity  gradient  darjdn  on  the  stator  surface,  which  is  obtained  by  calculating 
-  \l\}i{dplds)  in  Eq.  (5.1)  and  normalizing  it  by  where p  is  the  surface  pressure 

and  d^^  the  gap  between  the  rotor  and  stator  at  55  %R. 

The  contour  and  surface  plots  for  =  -10.30°  and  tz  =  10°  are  shown  in  Fig.  5- 
11(a).  From  Fig.  5-1 1(a),  it  is  apparent  that  the  active  region  is  concentrated  on  the  first 
10%C  of  the  stator’s  pressure  side.  Inside  this  region,  a  surface  vorticity  gradient 
d&Jdn  abruptly  rises  during  the  interaction  period,  as  we  observed  in  Fig.  5-5(d)  &  5- 
6(d).  Hence,  the  stator’s  loading  due  to  the  surface-vorticity  generation  is  highly  related 
to  the  noise  generation  sound  in  this  case.  It  turns  out  that,  during  the  interaction  period, 
a  strong  surface  vorticity  is  generated  into  the  flow  on  the  pressure  side  for  -22°  <  /?_^  <  0° 
at  six  stator’s  AoA. 

Since  the  surface-vorticity  generation  will  simultaneously  increase  the  stator’s 
lift,  a  transient  stator  response,  hence,  induced  by  the  vorticity  fluctuation  becomes  the 
source  to  radiate  the  sound  for  -22°  and  -15°  <  <  10°. 

5.4.3.  Cause  and  effect  of  rotor-stator  interaction  for  -22°  <  /?^  <  0° 

Therefore,  for  the  cases  that  -22°  <  <  0°,  the  cause  of  the  rotor-stator 

interaction  noise  is  as  follows:  During  the  interaction  period,  the  incoming  rotor’s  wake 
interacts  with  the  stator’s  boundary  layer  to  cause  a  sudden  pressure  drop  on  the  first 
10%C  of  the  stator’s  pressure  side.  And  hence  the  effect:  A  strong  vorticity  generation  by 
the  surface-pressure  drop  is  the  major  source  of  a  transient  stator  response  and  hence  the 
interaction-noise  radiation. 

A  schematic  description  of  the  noise-generation  mechanism  for  -22°  <P^<Q°  and 
-15°  <  a,  <  10°  is  illustrated  in  Fig.  5-12. 
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5.5.  The  Flow  Fields  for  Cases:  0°  <  J3^<  10°. 


In  the  rest  of  this  chapter,  we  are  about  to  survey  the  second  category  of  0°  < 
<10°.  The  =  9,61°  case  will  be  investigated  for  the  category  of  0°  <  /?^  <  10°. 

5.5.1.  Characteristics  of  surface  pressure  results  for  0”  <  <  10° 

In  Fig.  5-13,  surface  pressure  distributions  are  plotted  within  -60°  <  (jx  -12°,  for 
P^  =  9.61°  at  six  different  a/s.  We  find,  from  Fig.  5-13,  that,  at  each  a,  the  surface 
pressure  fluctuations  on  both  sides  of  stator  are  much  less  than  those  that  we  observed 
previously  in  the  p^  =  -10.30°  case  (-22°  <P^<  0°),  as  shown  in  Fig.  5-2.  Meanwhile,  the 
variations  of  surface-pressure  distributions  on  both  sides  for  -15°  <  a  <  10°  are 
apparently  too  small  to  distinguish  which  side  (or  both)  shall  be  the  active  region  for  the 
noise-generation  resource,  as  shown  in  Fig.  5- 13(a)  through  (f).  These  suggest  that, 
during  the  interaction  period,  the  noise  generation  may  not  be  consequent  on  the  stator’s 
surface-pressure  fluctuations  for  p^  =  9.61°  and  -15°  <  a  <  10°. 

In  other  words,  the  surface  pressure  variations  show  that  the  stator  response  is 
apparently  not  the  major  interaction-noise  resource  for  Q)°  <  p^<  10°  and  -15°  <  <  10°, 

which  is  exactly  disclosed  in  Chapter  4. 


5.5.2.  The  flow  fields  for  (f<p,<  10°  and  -15°  <  a.  <  10° 

The  velocity-contour  plots  for  p^  ~  9.61°  at  a,  =  10°,  0°  and  -10°  are  shown  in 
Fig.  5-14,  5-15  and  5-16,  respectively.  Generally  speaking,  since  the  rotor’s  effective 
AoA  p^  is  between  0°  and  10°,  the  rotor’s  wakes  observed  in  these  figures  are  apparently 
thinner  than  those  in  -22°  <  P^<  0°.  Moreover,  the  wakes’  angles  are  steep  and  close  to 
the  rotor’s  rotational  plane,  as  schematically  described  in  Fig.  5-17.  These  characteristics 
of  the  rotor’s  wakes  reveal  more  or  less  that  a  noise- generation  mechanism  for  the 
current  case  may  differ  from  that  in  -22°  <P^<  0°. 

An  overview  of  the  velocity-contour  plots  for  0°  <  >0,  <  10°  and  -15°  <  <  10° 

show  that  a  slight  variation  of  the  flow  field  in  the  vicinity  of  the  stator  surface  through 
the  entire  interaction  period  (Hu  and  Ho  1997).  This  agrees  to  the  surface  pressure  results 
in  this  category.  Since  the  noise  generation  still  occurs  during  the  interaction  period  from 
(t>  =  -60°  to  -30°  (for  example,  in  Fig.  3-2(c),  (d)  &  (e)  and  Fig.  3-3(e)),  this  implies  that 
the  interaction-noise  generation  may  be  from  the  a  transient  rotor  response. 
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5.6.  Summary 

The  physical  mechanism  of  the  interaction-noise  generation  for  -22°  <  /?  <  10° 
and  -15°  <  (S,<  10°  has  been  surveyed.  For  -22°  <  <  0°  and  -15°  <  <  10°,  the  flow 

interaction  is  based  on  the  first  10%C  of  the  stator’s  pressure  side  during  the  interaction 
period  from  ^  =  -60°  to  -30°.  This  interaction  causes  a  surface-vorticity  generation  from 
the  stator  into  the  flow;  hence,  a  transient  stator  response  is  induced  to  radiate  sound.  For 
0°  <  >9^  <  10°  and  -15°  <  >9  <  10°,  the  experimental  results  reveal  that  the  flow  interaction 
is  irrelevant  to  the  stator  response,  which  suggests  that  a  transient  rotor  response  may  be 
the  cause  of  interaction-noise  generation. 

As  a  summary,  if  y9^  decreases,  the  rotor’s  wake  will  be  tilted  more  toward  the 
stator  and  become  broader  in  width,  and  hence  the  strong  flow  interaction  will  take  place 
adjacent  to  the  stator  surface.  Namely,  for  a  fixed  a,,  the  stator  response  increases  as  /? 
-  increases,  and  it  decreases  as  decreases.  On  the  other  hand,  while  p^  decreases,  the 
majority  of  the  interaction-noise  generation  is  from  the  rotor  response,  as  discussed  in 
Chapter  4. 
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CHAPTER  6 
CONCLUSION 


The  noise  generated  from  a  rotor-stator  assembly  placed  in  an  incompressible 
pipe  flow  has  been  experimentally  investigated.  The  focus  of  this  study  is  to  understand 
the  physical  mechanism  of  the  interaction-noise  generation  by  correlating  the  acoustic 
field  to  the  unsteady  flow  field  around  the  rotor-stator  stage.  From  the  measured  noise 
signatures,  a  sudden  rise  in  the  noise  level  is  consistently  observed  and  identified  as  the 
interaction-noise  generation  from  -60°  to  -30°  in  phase  angle,  <j>,  for  -22°  <P^<  10°  and  - 
15°  <  a,  <  10°  (Chapter  3).  In  these  noise  signatures,  two  noise  patterns  can  be  easily 
distinguished  at  =  -10°  and  10°,  which  implies  that  there  exist  two  different 
mechanisms  for  the  noise  generation  in  the  survey  range. 

Based  on  the  Lighthill’s  acoustic  analogy.  Curie  (1955)  first  derived  the  sound 
field  radiated  by  a  stationary  solid  body  in  a  moving  medium.  Fujita  and  Kovasznay 
(1974)  simplified  Curie’s  result  to  study  the  rotor-stator  interaction  noise  generated  by  a 
transient  stator  response  in  a  free  space.  Since,  in  current  study,  the  fluid  flow  is  confined 
in  a  circular  pipe,  the  proportional  constant  in  Fujita-Kovasznay  formula,  Eq.(4.2),  is  no 
longer  valid.  Therefore,  we  can  only  calculate  the  source  term,  i.e.,  the  summation  of 


and 


R  dC^ 
dt 


in  Eq.(4.2),  from  the  unsteady  stator  response.  By  comparing  the 


calculated  source  term  with  the  measured  noise,  both  of  them  are  in  good  agreement  in 
peak  value  as  p^  decreases,  which  means  that  a  transient  stator  response  is  the  major 
resource  of  the  interaction  noise  generation  (Chapter  4).  On  the  other  hand,  as  p^ 
increases,  such  an  agreement  between  both  peak  values  decreases.  This  suggests  that  the 
unsteady  stator  response  contributes  less  to  the  noise  generation  as  p^  increases;  hence, 
most  of  the  noise  generation  are  from  a  transient  rotor  response.  In  other  words,  different 
mechanisms  for  the  interaction-noise  generation  results  from  a  variety  of  the  flow- 
incidence  angles  to  the  rotor,  i.e.,  P^.  This  is  because  p^  determines  the  flow  structure 
behind  the  rotor’s  wake  which  interacts  with  the  stator  to  generate  the  noise.  From  the 
velocity  results,  a  separated  wake  behind  the  rotor  observed  in  p^  <  0°  has  a  strong  flow 
interaction  with  the  boundary  layer  around  the  stator’s  leading  edge  during  the 
interaction-noise  generation.  At  this  interaction  period,  an  abrupt  decrease  of  surface 
pressure  on  the  first  10%C  of  the  stator’s  pressure  side  induces  unsteady  vorticity-flux 
production  into  the  fluid  flow,  which  causes  a  transient  stator  response  and  hence  sound 
radiation.  For  >  0°,  the  rotor’s  wake  is  thinner  and  tilted  to  the  rotor’s  side,  as 
observed  from  the  velocity-contour  plots,  and  the  experimental  results  show  that  small 
surface-pressure  and  velocity  fluctuations  occur  around  the  stator  during  the  interaction- 
noise  generation.  Accordingly,  the  transient  rotor  response  becomes  the  major  resource 
of  the  interaction-noise  generation  for  p^  >  0°. 
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Fig.  2-2  (  c  ) 


Fig.2-2  Arrangement  of  the  rotor  and  stator  (a),  definitions  of  phase  angle  ^  at  0°  and 
±180°  (b),  and  the  microphone  measurements  (c). 
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Rotor's 

Motion 


/?g  :  effective  AoA  of  rotor 
6  ;  angle  of  zero-lift  line 


Fig, 2-3.  Operating  condition  at  55%  radius  of  the  rotor  Note  that  both  and  are 

plotted  in  the  positive  sense. 
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condensor 


Fig.  2-4(a)  Configuration  of  a  transmission  line  to  the  stator  blade. 
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Plane  Waves 


Fig.  2-4(b)  A  conceptual  explanation  for  obtaining  the  frequency  responses  of 
transmission  line. 


^ref. 

Transmission 

^test 

Line 

Input 

Non-linear 

System 

Output 

Fig.2-4(c)  A  single-input  /  single-output  model  for  the  transmission-line  system. 


A 


test 


FIR  Filter 


Input 


Linear  System  Output 
( LTI  System). 


Fig.2-4(d)  Compensation  of  signal  by  FIR  filters. 


Fig.2-5  Configuration  of  the  hot-wire  probe  and  the  velocity  measurement:  a  typical 
measuring  region  on  the  55%R  cross-section  as  shown  in  Fig.  2-2(a). 
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Velocity  (  m/s  ) 


Frequency  (Hz) 

(b) 

Fig.2-7.  Gain  and  phase  responses  for  two  typical  pressure  holes:  diam.=  (a)  0.01 
(0.254mm)  and  (b)  0.02"  (0.508mm). 


0.010"  D  Pressure  Hole 


Frequency  (Hz) 
(a) 


0.020"  D  Pressure  Hole 


Frequency  (Hz) 

(b) 

Fig.2-8.  Compensation  of  gain  response  of  transmission  lines  by  FIR  filters  for  two 
typical  pressure  holes:  (a)  0.0  F’D,  (b)  0.02"D. 
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Sound  Pressure  (  Pa  ) 


Phase  Angle ,  4)  (  degree ) 


Fig.  3-l(c)  (  Uq  =  4.66  m/s  ,  55%R  =  12.23  m/s  :  pe=  5. 14°  ,  U(@  55%R  =  13.09  m/s.  ) 


-180  -90  0  90  180 

Phase  Angle ,  (>  ( degree) 


Fig.  3- 1(d)  (  Uq  =  4.66  m/s  ,  55%R  =  15.29  m/s  :  Pg=  9.05°  ,  U^  55%R  =  15.98  m/s.  ) 


1.20 


-180  -90  0  90  180 

Phase  Angle ,  ^  (  degree ) 

Fig.  3-l(e)  (  Uq  =  4.66  m/s  ,  55%R  =  19.1 1  m/s  :  pg=  12.30°  ,  U0  55%R  =  19.67  m/s.  ) 

Fig.  3-1.  Noise  signatures  measured  at  a  constant  freestream  speed,  4.66  m/s,  with  five 
different  rotor’s  speed  at  55%R:  (a)  6.1 1,  (b)  7.65,  (c)  12.23,  (d)  15.29  and  (e)  19.11  m/s, 
respectively.  The  resultant  local  incident  velocities,  f/, ,  are  (a)  7.68,  (b)  8.96,  (c)  13.09, 
(d)  15.98,  and  (e)  19.67  m/s,  while  their  corresponding  rotor’s  effective  AoA,  ,  are  (a) 
-11.33°,  (b)  -5.35°,(c)  5.14°,(d)  9.05°  and  (e)  12.30°,  respectively. 
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Fig.  3-2(e)  (  Uq  =  5.62  tn/s  ,  U,@  55%R  =19.11  m/s  :  Pe=  9.61°  ,  55%R  =  19.92  m/s.  ) 


Fig.  3-2.  Noise  signatures  measured  at  a  constant  freestream  speed,  5.62  m/s,  with  five 
different  rotor’s  speed  at  55%R:  (a)  6.11,  (b)  7.65,  (c)  12.23,  (d)  15.29  and  (e)  19.11  m/s, 
respectively.  The  resultant  local  incident  velocities,  V, ,  are  (a)  8.30,  (b)  9.49,  (c)  13.46, 
(d)  16.29,  and  (e)  19.92  m/s,  while  their  corresponding  rotor’s  effective  AoA,  ,  are  (a) 
-16.61°,  (b)  -10.30°,(c)  1.32°,(d)  5.82°  and  (e)  9.61°,  respectively. 
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Phase  Angle ,  <}>  ( degree ) 


Fig.  3-3(a)  (  Uq  =  8.43  m/s  ,  55%R  =  6. 1 1  m/s  :  Pe=  -28.07°  ,  55%R  =  10.41  m/s.  ) 


-180  -90  0  90  180 

Phase  Angle ,  <})  ( degree ) 


Fig.  3-3(b)  (  Uq  =  8.43  m/s  ,  Uj@  55%R  =  7.65  m/s  :  Pg=  -21.78°  ,  55%R  =  1 1.38  m/s.  ) 
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-3(c)  (  Uq  =  8.43  m/s  ,  55%R  =  12.23  m/s  :  Pe=  -8.58°  ,  55%oR  =  14.85  m/s.  ) 


-180  -90  0  90  180 

Phase  Angle ,  (()  ( degree) 

Fig.  3-3(d)  (  Uq  =  8.43  m/s  ,  U^  55%R  =  15.29  m/s  :  Pe=  -2.87°  ,  55%R  =  17.46  m/s.  ) 
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-180  -90  0  90  180 

Phase  Angle ,  ^  (  degree ) 


Fig.  3-3(e)  (  Uq  =  8.43  m/s  ,  55%R  =  19.1 1  m/s  :  Pe=  2.20°  ,  55%R  =  20.89  m/s.  ) 

Fig.  3-3.  Noise  signatures  measured  at  a  constant  freestream  speed,  8.43  m/s,  with  five 
different  rotor’s  speed  at  55%R:  (a)  6.11,  (b)  7.65,  (c)  12.23,  (d)  15.29  and  (e)  19.11  m/s, 
respectively.  The  resultant  local  incident  velocities,  U,  ,  are  (a)  10.41,  (b)  11.38,  (c) 
14.85,  (d)  17.46,  and  (e)  20,89  m/s,  while  their  corresponding  rotor’s  effective  AoA, 
are:  (a)  -28.07°,  (b)  -21.78°,(c)  -8.58°,(d)  -2.87°  and  (e)  2.20°,  respectively. 
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-180  -90  0  90  180 

Phase  Angle ,  (])  (  degree ) 

Fig.  3-4(a)  m/s  ,  55%R  =  6. 1 1  m/s  :  pe=  -35.47°  ,  U0  55%R  =  12.79  m/s. 


-180  -90  0  90  180 


Phase  Angle ,  (J)  ( degree ) 

Fig.  3-4(b)  {Vq^U.IA  m/s  ,  55%R  =  7.65  m/s  ;  Pe=  -29.76°  ,  55%R  =  13.60  m/s. 
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Phase  Angle ,  4>  (  degree ) 


Phase  Angle ,  ())  ( degree ) 

Fig.  3-4(d)  (  Uq=  11.24  m/s  ,  55%R  =  15.29  m/s  ;  Pe=  -10.32°  ,  Uf@  55%R  =  18.98  m/s.  ) 
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Fig.  3-4(e)  (  C/(j  =  1 1.24  m/s  ,  V,®  55%R  =  19. 1 1  m/s  :  Pe=  -4.46°  ,  f/^  55%R  =  22. 17  m/s. ) 


Fig.  3-4.  Noise  signatures  measured  at  a  constant  freestream  speed,  11.24  m/s,  with  five 
different  rotor’s  speed  at  55%R:  (a)  6.11,  (b)  7.65,  (c)  12.23,  (d)  15.29  and  (e)  19.11  m/s, 
respectively.  The  resultant  local  incident  velocities,  1),  ,  are  (a)  12.79,  (b)  13.60,  (c) 
16.61,  (d)  18.98,  and  (e)  22.17  m/s,  while  their  corresponding  rotor’s  effective  AoA, 
are:  (a)  -35.47°,  (b)  -29.76°,(c)  -16.58°,(d)  -10.32°  and  (e)  -4.46°,  respectively. 
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Freestream  Speed  ,  Uq  ( tn/s ) 


Fig.  3-5.  A  schematic  classification  by  on  a  plot;  along  Line  A,  =  -22°; 

along  Line  B,  =  0°;  and,  along  Line  C,  p^  -  10°.  The  cross  symbol  represents  the 
operating-condition  setup  for  noise  measurements,  while  the  circle  symbols  are  related  to 
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Fig.  3-6.  Peak  values  of  noise  v^.  the  resultant  speed,  t/,,  for  cases:  >  -22°.  For  a  fixed 
stator’s  AoA,  the  peak- value  data  are  curve-fitted  by  power,  log( Peak  Value)  =  b  \og(U) 
+  A  or  (Peak  Value)  =  af//,  where  power  b  represents  the  slope  of  each  line  in  the  log- 
log  plot.  Note  that  an  average  value  of  the  power  &’s  is  2  for  U,. 
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Uq  (  ni/sec  ) 

Fig.  3-7.  Non-dimensionalized  peak  values  of  noise  V5.  the  freestream  speed,  U^,  for 
cases:  >  -22°.  For  a  fixed  stator’s  AoA,  the  peak-value  data  are  curve-fitted  by  power, 
log(Peak  Value  /  q,)  =  d  \og(lJ„ )  C  or  (Peak  Value  /  q)  =  c  U/,  where  power  d 
represents  the  slope  of  each  line  in  the  log-log  plot,  where  q,  is  the  local  dynamic 

pressure  defined  as  ^  pU^  and  {/,  is  the  resultant  speed.  Note  that  an  average  value  of 

the  power  rf’s  is  1  for  U„. 
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Phase  Angle ,  0  ( degree ) 

Fig.  3-8.  Two  typical  noise  patterns  in  cases:  =  ±10°,  for  a  =  10°.  For  the  case  = 
-10°,  the  noise  reaches  a  valley  at  ^  =  -60°,  and  then  rises  abruptly  to  a  peak  at  ^  =  -30°. 
However,  for  the  case  (5^  =  10°,  the  noise  rises  gradually  from  a  valley  at  ^  =  -90°,  and 
then,  after  reaching  a  higher-noise  level  at  ^  =  -48°,  increases  rapidly  to  a  peak  at  ^  =  - 
30°. 
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Fig.3-9.  Noise  signatures  with  and  without  the  presence  of  stator  for  both  cases:  fie  =  (a) 
-10.30°  and  (b)  5.82°. 
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Fig.  4-2.  Time  traces  of  (a)  the  normal  force  coefficient,  Q,  and  (b)  the  tangential  force 
coefficient,  Q;  for  =  -21.78°  and  with  six  stator’s  AoA  a/s. 
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Fig.  4-3.  Time  traces  of  (a)  the  normal  force  coefficient,  Q,  and  (b)  the  tangential  force 
coefficient,  C^:  for  =  -10.30°  and  with  six  stator’s  AoA  a/s. 


Fig.  4-5.  Time  traces  of  (a)  the  normal  force  coefficient,  Q,  and  (b)  the  tangential  force 
coefficient,  C/.  for  =  2.20°  and  with  six  stator’s  AoA  a/s. 


75 


0.0015 


0.0010 


0.0005  -\ 

^  0.0000 
-0.0005 

-0.0010 

-0.0015 

-180  -90  0  90  180 

Phase  Angle ,  (j)  (  degree ) 


(b) 


Fig.  4-7.  Time  traces  of  (a)  the  normal  force  coefficient,  Q,  and  (b)  the  tangential  force 
coefficient,  Q:  for  =  9.61°  and  with  six  stator’s  AoA  a/s. 
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Fig.  4-14(c)  (0(3=  QO  ,  pg  @55%R  =  -21.78°  ,  55%R  =  11.38  m/s. ) 
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Fig.  4-14(e)  (  0(5=  -10®  ,  Pe  @55%R  =  -21.78®  ,  55%R  =  11.38  m/s. ) 


Fig.  4-14(f)  (  a5=  -15®  ,  p^  @55%R  =  -21.78®  ,  55%R  =  11.38  m/s. ) 

Fig.  4-14.  Measured  noise  V5.  source  term  :  for  =  -21.78°  with  a,  at  (a)  10°,  (b)  5' 
(c)  0°,  (d)  -5°,  (e)  -10°  and  (f)  -15°,  respectively. 
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Fig.  4-16(a)  (  0(5=  10®  ,  @55%R  =  1.32°  ,  55%R  =  13.46  m/s. ) 
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Fig.  4-17(a)  (  0(5=  10°  ,  @55%R  =  2.20°  ,  55%R  =  20.89  m/s. ) 
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93 


0.0075 


Measured  Noise 


0.0050 


0.0025 


I —  Source  Term 


S.P. 

1 — 771  0.0000 
jpf', 


-0.0025 


-0.0050 


-0.0075 


'  I  '  i  '  I  '  I  '  I  '  I  '  I  '  I  '■  I  '  I  '  I 

-180  -90  0  90  18 

Phase  Angle ,  (])  ( degree ) 

Fig.  4-17(c)  (  0(5=  0®  ,  pg  @55%R  =  2.20®  ,  Uf@  55%R  =  20.89  m/s. ) 
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Fig.  4-17(d)  (  0(5=  -5° ,  Pg  @55%R  =  2.20°  ,  55%R  =  20.89  m/s. ) 
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Fig.  4-18(c)  (a^=  0°  ,  @55%R  =  5.82°  ,  55%R  =  16.29  m/s. ) 
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Fig.  4-18(e)  (  0(5=  -IQO  ,  @55%R  =  5.82°  ,  55%R  =  16.29  m/s. ) 


Fig.  4-18(f)  ( cx5=  -15°  ,  Pg  @55%R  =  5.82°  ,  55%R  =  16.29  m/s. ) 

Fig.  4-18.  Measured  noise  v^.  source  term  :  for  =  5.82°  with  «  at  (a)  10°,  (b)  5°,  (c) 
0°,  (d)  -5°,  (e)  -10°  and  (f)  -15°,  respectively. 
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Fig.  4-19(a)  ( cxg=  10°  ,  pg  @55%R  =  9.61°  ,  55%R  =  19.92  m/s. ) 
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Fig.  4-19(b)  ( a5,=  5®  ,  pg  @55%R  =  9.61°  ,  55%R  =  19.92  m/s. ) 
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Fig.  4-19(f)  (  a  =  -15°  ,  @55%R  =  9.61°  ,  55%R  =  19.92  m/s. ) 

Fig.  4-19.  Measured  noise  vj.  source  term  :  for  =  9.61°  with  a,  at  (a)  10°,  (b)  5°,  (c) 
0°,  (d)  -5°,  (e)  -10°  and  (f)  -15°,  respectively. 


Fig.  4-21.  Amplitude  factor,  k,  vs.  interaction  angle,  y  :  plotted  for  two  constant 
freestream  speeds,  Ug ,  at  (a)  5.62  m/s  and  (b)  8.43  m/s,  respectively. 
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Fig.  5-l(c)  a  =-10°. 

Fig.  5-1.  4,500  hot-wire  stations  around  the  stator  from  the  leading  edge  to  32%  chord 
length  with  stator’s  AoA  a,  at  (a)  10°,  (b)  0°  and  (c)  -10°,  respectively.  The  stator’s 
chord  length  is  50  mm. 
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Fig.  5-2(e)  -6(P<(t)<-12°:  03= -10°,  pg= -10.30°,  9.49 m/s. 


Fig.  5-2(f)  -60°<<|)<-12°  :  as=-15°,  pe= -10.30°, 9.49 m/s. 

Fig.  5-2.  Non-dimensionalized  pressure  distributions  on  both  sides  of  the  stator  with 
stator’s  AoA,  a, :  (a)  10°,  (b)  5°,  (c)  0°,  (d)  -5°,  (e)  -10°,  and  -15°,  respectively,  for  = 
-10.30°  and  during  -60°  <  <  -12°.  Note  that  V,  =  9.49  m/s  with  f/,  =  5.62  m/s  and  t/, 

at  55%R  =  7.65  m/s. 
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Fig.  5-3.  History  of  velocity  contours  around  the  stator’s  leading  edge  during  -60°  <<f>^ 
-12°,  with  stator’s  AoA  a,  at  10°,  for  (5^  =  -10.30°  and  LJ,  =  9.49  m/s  at  55%R.  Note  that 
Vg  =  5.62  m/s  and  f/,  at  55%R  =  7.65  m/s. 
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Fig.  5-4(a)  Two  specific  lines  {Lines  A  and  B)  normal  to  the  stator  surface  at  the  leading 
edge  and  2%C  on  the  pressure  side,  respectively,  for  inspecting  the  stator’s  boundary- 
layer  variation. 
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Fig.  5-5  (a)  &  (b)  Temporal  variation  of  the  velocity  field  for  =  -10.30°  and  a,  =  10°; 
(a)  0.75  and  (b)  0.1  above  the  stator’s  leading  edge  along  Line  A  in  Fig.  5-4(a). 
Note  that  d  =  0.7 1  cm. 

sop 


Fig.  5-5  (c)  &  (d)  Time  trace  of  surface  pressure  and  vorticity  gradient  at  the  stator’s 
leading  edge  for  (i^  -  -10.30°  and  a,  =  10°. 
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Fig.  5-6  (a)  &  (b)  Temporal  variation  of  the  velocity  field  for  =  -10.30°  and  a,  =  10°: 
(a)  1.85  and  (b)  0.1  above  the  point  at  2%C  on  stator’s  pressure  side  along  Line  B 
in  Fig.5-4(a).  Note  that  d  =  0.71cm. 
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Fig.  5-6  (c)  &  (d)  Time  trace  of  surface  pressure  and  vorticity  gradient  at  the  point  of 
2%C  on  the  stator’s  pressure  side  for  =  -10.30°  and  a,  =  10°. 
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Fig.  5-7.  A  schematic  description  of  the  flow-interaction  process  for  /?,  =  -10.30°  and  a, 
=  70°:(a)a  narrow  flow  margin  formed  by  the  rotor’s  approach  in  -90°  <<f><  -48°,  and 
(  b  )  the  rotor’s  wake  impinging  on  the  stator’s  leading-edge  portion  during  -48°  <  ^  <  - 
30°. 
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Fig  5-8  History  of  velocity  contours  around  the  stator’s  leading  edge  during  -60°  <  ^  < 
-12^  With  stator’s  AoA  a.  at  0°,  for  p,  =  -10.30°  and  U,  =  9.49  m/s  at  55%R.  Note  that 
U,  =  5.62  m/s  and  V,  at  55%R  =  7.65  m/s. 
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Fig.  5-10.  A  local  coordinate  system,  s-n-z,  on  the  stator  surface. 
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Fig.  5-11  (f) 

Fig.  5-11.  The  contour  and  surface  plots  of  the  normalized  vorticity  gradient, 
^CO  *  ^  ^  5 

for  A  =  -10-30°  and  a  =  (a)  10°,  (b)  5°,  (c)  0°,  (d)  -5°,  (e)  -10°  and 

dn  U,  os 

(f)  -15°,  respectively.  The  dashed  contour  lines  show  the  values  which  are  less  than  zero. 
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Fig.  5-13(f)  -6CP<<t)<-12°  :  a^=-\5°,  (3g=  9.61°,  C/^  19.92 m/s. 

Fig.  5-13.  Non-dimensionalized  pressure  distributions  on  both  sides  of  the  stator  with 
stator’s  AoA,  a, :  (a)  10°,  (b)  5°,  (c)  0°,  (d)  -5°,  (e)  -10°,  and  -15°,  respectively,  for  p,  = 
9.61°  and  during  -60°  <<f)<  -12°.  Note  that  IJ,  =  19.92  m/s  with  t/„  =  5.62  m/s  and  [/,  at 
55%R=  19.11  m/s. 
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Fig.  5-14.  History  of  velocity  contours  around  the  stator’s  leading  edge  during  -60°  <  ^ 
<  -12°,  with  stator’s  AoA  a,  at  10°,  for  =  9.61°  and  U,  =  19.92  m/s  at  55%R.  Note 
that  U.  =  5.62  m/s  and  t/,  at  55%R  =  19.1 1  m/s. 
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Fig.  5-15.  (  Caption  on  next  page. ) 
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Fig.  5-15.  History  of  velocity  contours  around  the  stator’s  leading  edge  during  -60°  ^  <!> 
<  -12°,  with  stator’s  AoA  a,  at  0°,  for  =  9.61°  and  V,  =  19.92  m/s  at  55%R.  Note  that 
U,  =  5.62  m/s  and  C/,  at  55%R  =  19. 1 1  m/s. 
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Fig.  5-16.  ( Caption  on  next  page. ) 
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Fig.  5-16.  History  of  velocity  contours  around  the  stator’s  leading  edge  during  -60°  <  (j> 
<  -12°,  with  stator’s  AoA  at  -10°,  for  -  9.61°  and  U,  =  19.92  m/s  at  55%R.  Note 


that  U,  =  5.62  m/s  and  V,  at  55%R  =  19.1 1  m/s. 
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